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CALCULATED WIND-TUNNEL-BOUNDARY LIFT-INTERFERENCE 
FACTORS FOR RECTANGULAR PERFORATED TEST SECTIONS 

By Benferd L. Schilling and Ray H. Wright 
Langley Research Center 

SUMMARY 

Equations developed and presented in NASA Technical Report R-285 for approxi- 
mating the spanwise distribution of wind-tunnel -boundary interference on lift of wings in 
rectangular perforated-wall test sections are modified slightly to facilitate machine cal- 
culations and are then used to generate extensive tables of interference factors for a 
variety of wind-tunnel and wing parameters. Data are presented for horseshoe-vortex 
representations of small-span and large-span wings mounted at the center of rectangular 
test sections with five values of tunnel width-height ratio varying from 0.5 to 2.0 and with 
the ratio of permeability factor to compressibility factor ranging from 0.1 to 25.0. Span- 
wise distributions of upwash interference factor are given for each case. Use of the 
tables requires knowledge of the values of permeability factor applicable to the perforated 
test section for which lift interference factors are desired. Machine computer programs 
used in the calculations are presented as an appendix. 

INTRODUCTION 

During the last two decades, a number of wind tunnels have been constructed with 
test sections having perforated walls. Such walls are of particular advantage in testing 
through sonic speeds and at speeds slightly supersonic because they may reduce the 
severity of shock-wave disturbances reflected from the walls and impinging on the test 
models. However, these wind tunnels are commonly used also for subsonic testing, and 
if winged models are of appreciable size relative to the cross section of the tunnel at the 
test location, it is desirable to correct for the modification of test conditions due to the 
upwash interference of the tunnel boundaries. 

For constructional and operational convenience and for avoidance of focusing of 
reflected shocks in the transonic speed range, perforated-wall test sections are commonly 
made rectangular in cross section. An approximation method for estimating the wind- 
tunnel -boundary upwash interference along the span of a lifting wing mounted at the center 
of such a perforated-wall test section was developed in reference 1. However, the exten- 
sive calculations required to obtain numerical values of upwash interference factors 



involve the numerical evaluation of multiple infinite integrals containing integrands 
singular at the zero point of the domain of integration. In order to facilitate the applica- 
tion of the theory of reference 1, machine computing programs for calculation of upwash 
factors were prepared and are presented herein along with tables of upwash factors 
applicable to a range of practical configurational and operational parameters. Figures 
constructed from the tables illustrate the variation of upwash interference factor with 
effective permeability of the test-section walls, with spanwise location along the wing, 
and with ratio of width to height of the test section. Although the theory of reference 1 
provides for the possibility that the permeability of the top and bottom walls is different 
from that of the side walls, it is assumed uniform for these calculations. 

SYMBOLS 


A cross-sectional area of test section 

a variable used for convenience to represent more complicated expression in 

equation (A12) 

b semiwidth of test section 


0^ lift coefficient 


Fi(q),F{(q) 

G(q,r),G’(q,r) 
Gi(q,p),G-(q,r),G'i (q,r) 
Ki(q,r),Ki(q,r) 
H(q,r,p),H'(q,r,p) 
Hi(q,r,p),Hi(q,r,p),H'i'(r,p,0) 


\ 


Symbols used for convenience to represent 
functions of the indicated dummy variables 


Hy(q,r,p),Hy(q,r,p),Hy(r,p,0) 
Li(q,r,p),L{(q,r,p) 
Ly(r,p),Ly(r,p) ) 


h semiheight of test section 


M Mach number 


2 


p,p' dummy variables of integration 

q dummy variable of integration 

R permeability factor 

r dummy variable of integration 

5 area on which lift coefficient is based 

s semispan of horseshoe vortex representing wing 

V velocity of tunnel test stream 

v upwash interference velocity, positive in direction of Z -axis 

X,Y,Z axes of rectangular Cartesian coordinates 

x,y,z rectangular Cartesian coordinates, x in direction of tunnel flow, y along 

direction of wing span, and z vertical 

0 = \jl - M 2 
r circulation 

6 upwash interference factor 

6 dummy variable of integration 

p dummy variable of integration 

Subscripts: 

The subscript j is used as a secondary subscript only as a bookkeeping device. 
The subscript i is used as a primary subscript with the following definitions: 

2 pertaining to vertical boundaries 

3 pertaining to horizontal boundaries 
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4 pertaining to effect of horizontal boundaries on interference potential inside 

test section due to vertical boundaries 

5 pertaining to effect of horizontal boundaries on interference potential outside 

test section due to horizontal boundaries 

UPWASH INTERFERENCE APPROXIMATION EQUATIONS 

Approximation equations for calculating the upwash interference velocity v due to 
the boundaries of a perforated wind tunnel along a wing mounted at the center were derived 
in reference 1. The wind tunnel was assumed rectangular with semiheight h and semi- 
width b, and the wing was represented by a horseshoe vortex with span 2 s and circula- 
tion T. A rectangular Cartesian coordinate system was used with coordinates z posi- 
tive in the direction of lift, y lying along the wing span, and x positive downstream in 
the direction of the trailing vortices. A schematic diagram showing the relationships 
between the various parameters is given in figure 1 . 

The equations of interest in reference 1 are equations (A14), which gives the inter- 
ference velocity due to infinite vertical boundaries; equation (BIO), which gives the inter- 
ference velocity due to infinite horizontal boundaries; and equations (CIO) and (C12), 
which give increments of interference velocity to satisfy more nearly the boundary condi- 
tions when both horizontal and vertical boundaries are present. These equations involve 
the permeability factors R 2 along the vertical boundaries and R 3 along the horizontal 
boundaries. If R 2 = R 3 = R; that is, if the permeability factor is assumed uniform and 
equal on all boundaries, then the equations approximating the various components of the 
interference velocity can be written in terms of R//3 (where /3 = ^1 - m 2 and M is 
Mach number) instead of separately in terms of R 2 , R 3 , and /3. 

The upwash interference factors 6 j, which are obtained from their respective 
upwash interference velocities by the relationship 

1 SVCl 


can be put in the form 


«i 


b/h h 

s7h T 


from which the total upwash interference factor is calculated as 

i 


( 1 ) 

(2) 
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The four equations referred to in reference 1 have been rewritten by means of 
equation (1) to give their respective upwash interference factors. These have also been 
rewritten so that the terms R and |3 never appear singly, but always appear in the 
form R//3. Finally, the equations have been rearranged and transformed to the extent 
necessary to present them in forms more readily adaptable to machine calculations. The 
details of these manipulations are given in appendix A. The machine computer programs 
are presented in appendix B. 


RESULTS AND DISCUSSION 

As mentioned in the preceding section, the approximation equations have been 
rewritten so that the upwash interference factors are expressed in terms of R//3 rather 
than singly in either the permeability factor R or in /3. For the purposes of calcula- 
tion, j3 was allowed to take the values 1.0, 0.8, 0.6, 0.45, and 0.3 for each of the values 
of R, which were 0.1, 0.45, 2.0, and 7.5. Thus, upwash interference factors were calcu- 
lated for 20 values of R//3 ranging from 0.1 to 25.0. The calculations were made for 
this range of R//3 for each of five different values of tunnel width-height ratio b/h, 
which were 0.5, 0.75, 1.0, 1.5, and 2.0, for both a small-span wing (s/b = 0.3), where s/b 
is the ratio of wing span to tunnel width, and a large-span wing (s/b = 0.7). Finally, the 
interference factors were calculated for three locations along the wing span (y/s = 0.0, 

0.5, and 1.0) for each value of the parameters s/b, b/h, and R//3. 

The resulting calculated upwash interference factors for a small-span wing mounted 
in the center of a rectangular perforated wind tunnel are presented in table I. Each page 
of the table presents the variation of upwash interference factors along the wing span for 
the entire range of R//3 for one of the five values used for b/h. Table II presents the 
corresponding calculations for a large-span wing. Values are given for each of the indi- 
vidual upwash interference factors 63, 63, 64, and 65 (as defined by eq. (1) in con- 

junction with subscript definitions in the list of symbols) as well as for the total upwash 
interference factor 6. 

There is always the possibility of calculation errors (such as round-off errors) in 
calculations of this type, and a certain amount of effort is required to ascertain the extent 
of these errors in order to decide what level of confidence to place in the final calcula- 
tions. The infinite integration limits in these particular calculations necessitated an 
arbitrary truncation of the calculating procedure beyond a certain point in order to con- 
serve computer time. To some extent with the double integrals, and especially with the 
triple integrals, compromises had to be made between the level of accuracy desired and 
the calculating time used. Therefore, the integration range was broken up into small 
equal intervals the lengths of which were chosen from information obtained in preliminary 
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computer runs of the program. Integrations were then made successively over each of 
the intervals, the integration limits increasing by the interval length with each integration 
step. When the integral over the last interval was less than some arbitrarily chosen 
small number (a change of one unit in the fifth decimal place in this case), the integral 
was assumed to have converged and the integration procedure was terminated. In these 
calculations, the rapid approach of the integrands toward zero with increase of the inte- 
gration variables lent credence to this criterion of convergence. 

The accuracy of machine calculations also depends upon the number of points per 
integration interval used in the machine integration routine (the Gaussian quadrature 
method was used for these calculations). Here, again, compromises had to be made 
between desired accuracy and calculating time. For the calculations herein, the number 
of points used was increased over several preliminary runs of the program, and the final 
number of integration points chosen was that for which the value of the integral changed 
by no more than one unit in the fifth decimal place between successive trials with - 
increasing numbers of points. 

Special consideration and handling must be given to integrals whose integrands con- 
tain terms causing the integrands to oscillate about zero as the integration variable 
increases. In several instances in the integrals reported herein (as discussed in 
appendix A), the integrands contained products of sines and cosines whose arguments 
contained the same variable but differed in the constants by which the variable was multi- 
plied. Thus, several terms in the integrand were oscillating about zero at different fre- 
quencies. This behavior indicates the possibility that oscillations might combine within 
some particular integration interval in such a way as to cause the integral over that inter- 
val to be very small. The calculation process might then be truncated prematurely 
because of a spurious indication of convergence of the integral. The approach used in 
handling this problem was to combine the sine-cosine product terms by means of multiple- 
angle trigonometric identities into sums of sine and eosine terms whose arguments were 
the same and then to choose the integration interval to extend over some multiple of a 
half-cycle of the new argument. Tables I and n present data only to three decimal places, 
but the tests on convergence and number of points per integration interval were to five 
decimal places; therefore, the values in the tables are assumed free of arithmetic trunca- 
tion errors. 

If the approximation equations for the various components of upwash interference 
factors are valid, then the total upwash interference factors should approach the values 
for the completely closed tunnel as R//3 approaches zero and the value for the com- 
pletely open tunnel as R//3 becomes very large. Such behavior is confirmed by the 
curves of figure 2 which show the variation of total upwash interference factor at the 
center of a small-span wing mounted in the center of a rectangular perforated wind tunnel 
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as a function of R/j3 for several values of tunnel width -height ratio b/h. Figure 3 
shows the corresponding curves for a large span wing, and again the calculated values 
tend toward the closed-tunnel values for small R//3 and toward the open-tunnel values 
for large R//3. Note that on the logarithmic scale of figures 2 and 3, the point on the 
abscissa corresponding to the closed tunnel (R/fi = 0) lies infinitely far to the left, and 
the point corresponding to the open tunnel (R//3 — °°) lies infinitely far to the right. 

The open-tunnel and closed-tunnel upwash interference factors at the center of the 
test section, which are given for comparison, were calculated from equations independ- 
ently derived by the method of images as in reference 2. For the closed test section, 6 
at the center of the wing is 



and for the open test section, 6 at the center of the wing is 



Machine computer programs corresponding to these equations are given in appendix B. 

Plots and cross plots made from the extensive data presented in tables I and II 
afford a great deal of insight into the variation of upwash interference factors with vari- 
ous tunnel and wing parameters. No attempt will be made in this paper to present a com- 
prehensive set of figures showing the interrelationships of the upwash interference fac- 
tors with the various other parameters. Two figures (figs. 4 and 5) illustrative of the 
type of information available are given. Figure 4 shows the variation of total upwash 
interference factor along the spans of both small-span and large-span wings for several 
values of R//3 and for several tunnel width -height ratios. Figure 5 shows the variation 
of upwash interference factor with tunnel width-height ratio for several values of R//3 
for both the small -span and large -span wings. 

A comparison of figures 4(a) and 4(b) shows that the variation of the interference 
along the span is small for the small -span wing but substantial for the wing spanning 0.7 
of the test-section width. Also, the largest spanwise variation occurs for the test sec- 
tion having the smallest width-height ratio, an effect that might have been anticipated 
since in this case a large part of the interference is produced by the side walls. Note 
however, that in all cases for which the permeability is such as to produce nearly zero 
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interference at the center of the wing, the interference at the wing tips is also nearly 
zero. 

Figure 5 shows that except for the nearly open tunnel with height greater than the 
width, the downwash interference at the center of a wing spanning a given fraction of the 
width of the test section increases with increase of the test-section width-height ratio. 

CONCLUDING REMARKS 

Equations developed in NASA Technical Report R-285 for approximating the span- 
wise distribution of wind-tunnel -boundary interference on lift of wings in rectangular 
perforated-wall test sections have been modified to facilitate machine calculations and 
then used to generate tables of interference factors for a variety of wind-tunnel and wing 
parameters. 0 

The upwash interference factors presented in the tables are applicable only if the 
permeability factor R is known. Inasmuch as R depends not only on the geometry of 
the perforated walls but also on operating conditions such as Reynolds number, Mach 
number, and boundary-layer thickness, it cannot be calculated but must be experimentally 
determined. The effective permeability factor may also vary from place to place over the 
perforated walls so that the use of some average values may be required. 

The approximation method used in the calculations is believed to yield upwash inter- 
ference factors of adequate accuracy. Because the wing is represented by a single horse- 
shoe vortex, the calculated upwash interference factors are applicable only to wings of 
reasonably small chord. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 16, 1969. 


L 

£ 

! 


8 



APPENDIX A 


ADAPTATION OF UPWASH INTERFERENCE EQUATIONS 
TO MACHINE CALCULATIONS 


Equation (A14) of reference 1, which represents the upwash interference velocity 
due to infinite vertical boundaries, can be rewritten by means of equation (1) in the body 
of the present report to give the upwash interference factor 


where 


and 


6 2= JLb^l 
2 2 it 2 s / h 


\ F2(q)dq - - \ \ G2(q,p)dp dq 

JO R Jq Jo 


F 2 (q) = 


7T sinh^q^cosh^q^ 


h q 

e h sinhl^q 


(Al) 


G2(q,p) = 


q 2 sinh(^/3 2 p 2 + q 2 jcoshf^/3 2 p 2 + q 2 


//3 2 p 2 + q 2 


/3 2 p 2 + q 2 s inh2|b^2p2 + q2 j + Tp2cosb?(^Jffi]}2 + q2 j" 


If the transformation of variables 


jSp = r cos 6 


q = r sin 6 


is made in the double-integral portion of equation (Al), then the differential element 
/3 dp dq is replaced by r dr d0 and the double integral becomes 


M 

R 


’ sinh^ r^ cosh^ r^ 


0 


>tt/2 


sin z 0 d0 dr 


cosh 2 ^) J 0 (|) 2 tanh 2 ^r) + cos 2 6 


Rearrangement of the integral over 9 by means of the double -angle trigonometric 
identities gives an integral in 9 which can be integrated analytically. If the dummy 
variable of integration r is replaced by q, the result can be combined with the single - 
integral portion of equation (Al) to give the upwash interference factor 
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52 = SS$I 0 " F2(q)d ' 1 (A2) 

where 



In the limits, as the permeability factor R approaches zero (closed-tunnel case) and 
infinity (open-tunnel case), equation (A2) approaches the same limits as those obtained in 
reference 1. 

The upwash interference velocity due to infinite horizontal boundaries, given by 
equation (BIO) of reference 1, is written in terms of the upwash interference factor as 


where 


* 1 b/h 

53= 2^S7h 


n°° rt rt°° 

J o F 3 (q)dq - - J G 3 (q,p)dp dq 


F 3 (q) = 


7T cos^q)sin(|qj 


cosh q 


(A3) 


and 


cos 


G 3 (q,p) = 




^2 cosh 2 (\Jii 2 p 2 + q 2 j ■ 


/3 2 p 2 + q2 


sinh 2 ^/3 2 p 2 + q2 


In order to write the second integral in terms of R//3, the transformation p = r/R 
(where R is a finite constant) is made. The machine calculations are made more 
accurate if the product cos^q)sin^q) is written as 



Finally, using the identity sinh 2 x = cosh 2 x - 1 and simplifying give equation (A3) in the 
form 
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« 1 bA 

63 = i^S7h 


r oo r oo /-»oo 

F3(q)dq - ^ ^ G3(q,r)dq dr 


(A4) 


where 


and 


f i(q> - f 


» sln [(f + h) q 1 tsln [(S'h) q 


e^ cosh q 


G3(q,r) = 


(fr ♦ q2 ){; ln [(g + h) q ] * sln [(f ' 


^ + ^ - q2 c ° sh2 (/^ + q2 ) _ 4 


The upwash interference factors due to the effect of horizontal boundaries on the 
interference velocity potential inside the test section due to vertical boundaries 64 
(from eq. (CIO) of ref. 1) and to the effect of horizontal boundaries on the interference 
velocity potential outside the test section due to horizontal boundaries 65 (from 
eq. (C12) of ref. 1) are given by 


k _ 1 b/h 
4 " tt 2 s/h 


-C (* G(q,r)G 4 (q,r)dq dr + f C f H(q,r,p)H 4 (q,r,p)dq dr dp 

J 0 Jo Jq j 0 j 0 


(A5) 


and 


«. 1 b/h 

6 5“^2i7h 


r»oopoo rt n°° r*°° r*°° 

\ \ G(q,r)G5(q,r)dq dr + - 4 ~ \ \ \ H(q,r,p)Hs(q,r,p)dq dr dp 

Jq j 0 j 0 j 0 j 0 


(A6) 


where the symbols common to both equation (A5) and equation (A6), that is, G(q,r) and 
H(q,r,p) represent 


G(q,r) = 


sinh^qjcos^rjcos q 




(q2 + r2^e^ cosh r 
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and 


H(q,r,p) = 


q sinh^ \Jfi 2 p 2 + q 2 jcos^r) 


(gV ♦ + r2 + i cosh2^V7^e^^ ? ) 


The symbols G4(q,r) and H4(q,r,p) in equation (A5) represent 


Gi(q,r) = 


q sinh^qjcos^rj + r cosh^q)sin(^r) 
stahgq) 


and 


H4(q,r,p) = H 4)1 (q,r,p)jH4 > 2(q,r,p)H 4j 3(q,r,p) - H 4j4 (q,r,p)H4 5 5(q,r,p)] 


where 


H 4,l(q*r,P) = 


^3 2 p 2 + q 2 sinh^y^ 2 p 2 + q 2 jcos(^rj + r cosh^^p 2 + q 2 jsin^ 


sinh 2 ^ys 2 p 2 + q 2 j+ p 2 co6h 2 ^y8 2 p 2_ + _ q^ 


H4,2(<l> r >P) = 


sinh^y/3 2 p^ + q 2 j + coshf^y^p 2 + q 2 j 


//3 2 p 2 + q 2 


p 2 sin q sinh//3 2 p 2 + r 2 Q « 

H 4 ? 3 (q,r,p) = + ^ cos q cosh + r* 

y/3 2 p 2 + r 2 


sinh(^i//3 2 p 2 + q 2 ) p 2 coshfe\//3 2 p 2 + q 2 

H 4 4 (q,r,p) = — ^ 1 - 

R 2 £ 2 p 2 + q 2 


and 


, . , nr? — ? q cos q sinh i/s 2 p 2 + 1-2 

H4 > 5(q,r,p) = sin q cosh + r z == 

03 2 p 2 + r 2 
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In equation (A 6 ), the symbols G 5 (q,r) and H 5 (q,r,p) represent 

Gfj(q,r) = q cos^r ) - r sin(|r ) 


and 


H 5 (q,r,p) 


H 45 (q,r,p) 

j3 2 p2 + q2 



+ 




Consider for a moment only the double -integral portions of equations (A5) and (A 6 ). 
These represent the case of the completely closed tunnel. As with equation (A3), the cal- 
culation procedures are made more accurate by writing the product cosQ^rjsin^r^ and 

the product cos^rjcos^rj as sines and cosines with arguments (j" + ^j r an d _ 

The double integrals of equations (A5) and (A 6 ) can thus be written 


and 


where 


64 (double -integral portion) = — ^ ^4^ C C G'(q,r)G 4 (q,r)dq dr 

2 7 r2 s /h ^0*^0 


65 (double -integral portion) 


G ’ (< i" r > G 5«i> r ) d <‘ dr 


G'(q,r) 


G 4 (q,r) 


K 4 (q,r) 

K4(q,r) 


cos q sinh^q 

(q2 + r 2j e h COS h r 

K 4 (q,r) + K^r) 


slnhgq) 




= q S inh(|q)cos 

(H)3 

+ r cosh^q^sin 

li - Di 

= q sinh^q^cos 

(HDi 

+ r cosh ^qj sin 



(A7) 


(A 8 ) 
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G’ 5 ’(q,r) = K5(q,r) - K 5 (q,r) 


4(q,r) = q cosgg - i)rj - r S ln{(£ - 2)r] 

Kg(q,r) = q «*[(£ + £)r] - r sin |(g + *)rj 

The triple -integral portions of equations (A5) and (A6) are formulated in terms of 
R//3 by writing the dummy variable p as p'/R. If the obvious cancellations are made 
and the expression for H4 j2 (q,r,p) is simplified by writing it in terms of exponentials 
rather than as the sum of hyperbolic terms, triple integrals in terms of the dummy vari- 
ables q, r, and p' are obtained. For convenience, the prime on the variable p* is 
dropped, and the triple -integral portions of equations (A5) and (A6) become 


and 


64 (triple-integral portion) = f f f H' (q,r,p)H4(q,r,p)dq dr dp (A9) 

7 t j s/h Jq Jo j 0 

65 (triple -integral portion) = J H* (q,r,p)H5(q,r,p)dq dr dp (A10) 


where 


H’(q,r,p) = 


q cos(gr)sinh^^p7 ^ 


+ q 2 


W R2 W + q2 + r 2) 

R 2 ) 


Lfi sinh 2 + r 2 


\ R 2 


£si + r 2 


+ cosh 2 /^ + r 2 
V R 2 


/ 


H4(q,r,p) = H4 >1 (q,r,p)jH4 j2 (q,r,p) - H4 >3 (q,r,p)H4 j4 (q,r,p) 


H4 ?1 (q,r,p) = 


^ + q2 + r cos1, (ki 

(r^ + q2 ) s “ 2 (ki/^t + q2 ) * p2cosh2 |/^ + q2 ) 
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H4 j2 (q,r,p) = 


+ r2 i gg-^ 

r + q cos q cosh p- g - + 


fS ^ 2 


^V + r 2 

R2 


r2 


H 4 ? 3 (q,r,p) = sinh^J^|^ + q 2 


p 2 cosh(£ r 2 P 2 + q2 

Vfty 

g 2 p 2 o 

li- + q 


H 4 )4 (q,r,p) = sin q coshj£-|- + r 2 - 


S3 ; q cos q sinh l&j* + r 2 


« + r 2 

R 2 


and 


. H 4 4 (q,r,p) 
H 5 (q,r,p) = — 2 


ff 2 p 2 

R 2 


+ q* L7 



cos(|r) - r Bin(gr) 


Investigation of the integrand of equation (A9) reveals a singularity at the simulta- 
neous zero of all three integration variables. For small values of R//3 this singularity 
can be ignored in the numerical integrations, but for moderate and large values it leads 
to difficulties in the calculation. In order to carry out the machine calculations for the 
entire range of R//3 considered herein, it is necessary to break the triple integral of 
equation (A9) into several integrals, one of which is used to obtain an approximation for 
the integral near the simultaneous zero of the integration variables, and the others are 
used to calculate the integral over the remainder of the integration range. Thus in equa- 
tion (A9), assume that the integration variables p, q, and r are all so close to zero 
that the transcendental functions in the integrand can be well approximated by the first- 
order terms in their series representations. After the indicated substitutions and sim- 
plifications are made, equation (A9) becomes 
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64 (triple -integral portion) 


p— 0 

q— 0 
r— 0 






q^dr dq dp 




>e 2 P € 3 q 2 dq dp 

'0 ^0 fe + q2 ) + (^h) 2 


(All) 


where the integration limits ej_, €3, and 63 must be chosen small enough to make the 
approximation to equation (A9) valid. Equation (All) becomes more amenable to integra- 
tion by use of the transformation of variables 


^ = p sin 6 


q = p cos 6 


Whence, 




R 

and the differential element is — p dp 68. As a result of applying this transformation 

/3 

and using the identity sin 2 0 + cos 2 0 = 1, equation (All) becomes 


64 (triple -integral portion) 


near 

origin 


= £fl 5 C* 2 p dp C 2 cos 2 6> 68 
3 P ^0 ^0 a 2 sin^0 + p 2 cos 2 0 


(A12) 


where a 2 = p 2 + . In this expression, the integration limit . ejj is the equivalent 

in the transformed coordinate system of €3 and €3 in the former system, but the inte- 
gration region is now a quarter of a circular disk instead of a rectangle. Equation (A12) 
is a standard integral form in 8 which can be integrated to give 


64 (triple -integral portion) 


near 

origin 
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This expression may then be integrated over p to yield 


64 (triple -integral portion) 


near 

origin 


e l (b/h) 2 
2tt 2 R/0 



+ 



(A 13 ) 


In order to perform the integration of equation (A 9 ) over the remainder of its inte- 
gration range, it is necessary first to apply the same coordinate transformation used to 
obtain equation (A 12 ). The following equation results: 


where 


64 (triple -integral portion) = [Eq. (A 13 )| + -tt ^ 4 ^ C dr C dp C H4(r,p,6»)d0 

jro s/n Jq j 6 > Jq 




p°0 r£ ^ P2 

+ \ dr \ dp \ HV(r,p,0)d0 

J e 1 J 0 J 0 


E 
2 


Hj(r,p, 9 ) = Hi' ;1 (r,p,e)H4 j2 (r,p,9)[H4 ; 3(r,p,e) + H4 ;4 (p,0)H4 ;5 (r,p,<))j 

!inh (pP) 


cos 6 sinh 


H4,l(r,P,») - b 


^ ^p2 + r 2 ^ sinh 2 ^pj + sin 2 0 cosh 2 ^p^ 


H 4 j2 ( r ,P, 6 ) = 


L4 )2 ( r >P) + ^4,2 ( r >P) 


/R\ 2 p 2 si n 2 e sinh 2 /p 2 sin 2 6 + r 2 + cosh 2 ^ p 2 s i 
V* / p 2 sin 2 0 + r 2 


sin 2 0 + r 2 


(A 14 ) 


L^ ;2 (r,p) = p slnh( 5 p)co S |(£ - *)rj + r co S h(£p) S m|(£ - 
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/b„\ „ 

/b y\ 

(k p ) cos 

|ve + J) g 


+ r cosh 


(ii p ) sin (k + ?) 


H4 5 3(r,p,e) = e 


-P 


R \ 2 p sin 2 0 sin(p cos 0)sinh\/p2sin20 + r 2 

I. 


yp2sin^0 + r 2 

+ cos 6 cos (p cos 0)cosh^p 2 sin 2 0 + r 2 


H 4 )4 (p, 0 ) = sinh^pj sin 2 0 cosh^pj 


and 


TT t» , p cos 0 cos(p cos 0)sinhl/p 2 sin 2 0 + r 2 

H 4, 5 fr,P,e) ^ 

)jp* sin* 6 + 


sin(p cos 0)cosh^/p 2 sin 2 0 + r 2 


In the expression for H 4> 2(r,p,0), the terms and L 4j 2 (r,p) are the result of 

the product cos^r^sin^r^ and the product cos ^rjcos^rj, similar to the analagous 
expressions of equations (A7) and (A8). In a like manner, the integrand of equation (A10) 
may be written as 


H T (q,r,p)H 5 (q,r,p) = H 5jl (q,r,p)H 5 }2 (q,r,p) 


where 

H 5,l(q»r,P) 
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Hg,a(<l>r.P) = H i’ 4 2 (C1 ’ r,r - [L5(q, r ,P) + I- 5 (q,r,p) 


Ls^^jP) 

and 

L5<q,r,p) 




APPENDIX B 


FORTRAN PROGRAM FOR CALCULATING SPANWISE VARIATIONS IN 
W I NO-TUNNEL- 8 CUNC ARY LIFT-INTERFERENCE FACTORS FOR WINGS OF 
VARYING SPAM CEN TER-MCUNTEC IN RECTANGULAR PERFORATED TEST 
SECTIONS OF VARYING W IDT H-TC- HEIGHT RATIOS* 


THIS PROGRAM HAS WRITTEN IN CDC FORTRAN , VERSION 2.1, TO RUN ON CCC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY CN TEE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY 0 


IN ORDER TC CONSERVE CCNFLTER TIME AND SPACE, THIS PROGRAM 
IS WRITTEN AS A NUMBER OF SMALLER, INDEPENDENT PROGRAMS WHOSE 
PUNCH! D-CARD OUTPUTS ARE COLLATED INTC FINAL FORM FOR PRODUCTION 
OF TABLcS BY MEANS OF APPROPRIATE COMPUTER COLLATING PROGRAMS. 

FOR THE SAKE CF CONSISTENCY, TERMS COMMON TO ALL PROGRAMS HAVE 
BEEN ASSIGNED THE SAME NAMES IN ALL PROGRAMS. THESE TERMS FOLLOW 


A P , AC, AR = LOWER LIMITS CF INTEGRATION ON THE INTEGRATION 
INTERVALS OVER The DUMMY VARIABLES P, C, AND R . 

BH = A CNE-DI ^ENSICNAL ARRAY IN WHICH ARE STORED THE VALUES CF 
TUNNEL SEMI — WICTH TO SEMI-HEIGHT RATIO, ' B/H 

BP, BG, BR = UPPER LIMITS CF INTEGRATION CN THE INTEGRATION 
INTERVALS OVER THE DUMMY VARIABLES P, C, AND R. 

8 S THE RECIPROCAL OF SB. 

ST = A ONfr-DI MENSICNAL ARRAY IN whICH ARE STORED THF VALUES OF 
BETA TO BE USED IN THE C A LC 0 L A T I C N S • 

CNV = THE NUMBER WHICH DETERMINES INTEGRAL CONVERGENCE. IF THE 
EVALUATION OF THE INTFGRAL CVER A PARTICULAR INTERVAL IS LESS 
THAN CNV, CONVERGENCE IS ASSUMFD ANC THE INTEGRATION PROCESS 
TERMINATED. 

FP, FQ, FR = ONE DIMENSIONAL ARRAYS IN WHICH ARE STORED THE VALUES 
CF INTEGRAND EVALUATIONS CURING THE INTEGRATION PROCESSES. 

FUNCP, FUNCQ , FUNCR = NAMES CF SUBROUTINE SUBPROGRAMS WRITTEN TO 
EVALUATE THE INTEGRANDS IN THE DUMMY VARIABLES P, Q, AND R • 

IB OR THH = INDEX CN THE VARIABLE BH 

INFINTP, INF INTO, INFINTR = NAMES CF SUBROUTINES WHICH USE MGAUSP, 
MGAUSQ, AND MG A J SR » RESPECTIVELY, TO EVALUATE THE ^INFINITE* 
INTEGRALS OVER THE DUMMY VARIABLES P, Q, AND R. 

IP = INDEX ON VARIABLE RP . 
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IS 

Ok 

I Sh = 

INDEX 

ON 

THE 

VARIABLE 

SH 

I Y 

OR 

IYH = 

INDEX 

ON 

THE 

VARIABLE 

YH 


JP, JO* JR = THE MAXIMUM NUMBER OF INTEGRATION INTERVALS USED BEFORE 
NCN-CONVEPGENCE CF THE INTEGRALS IS ASSUMED. 

MGAUSP, MGAUSC, MGAUSR = NAMES CF SUBROUTINES FOR EVALUATING 
INTEGRALS BY THE GAUSS QUADRATURE METHCC CVER THE DUMMY VARIABLES 
P, C, AND R * RESPECTIVELY. 

NFP, NFC, NFR = INTEGERS DEFINING THE NUMBER OF INTEGRANDS TO BE 
EVALUATED IN THE INTEGRATION SUBROUTINE CALL FOR ANY CF THE 
DUMMY VARIABLES P, 0, OR R. 

NMY = NUMBER CF POINTS ALONG THE WING SPAN AT WHICH CALCULATIONS 
ARE TO BE MADE 

NP, NO * NR = INTEGERS DETERMINING THE NUMBER OF POINTS THE GAUSS 
CUADRATURE PRCCECURE USES PER INTEGRATION INTERVAL. THE NUMBER 
OF POINTS USED IS TEN TIMES THE VALUE CF NP, NC, OR NP. 

PI = 2.14154 it 5 26' 

PNC, CMC , RNC = VARIABLES WHICH SET THE INTERVAL LENGTH FCR 
INTEGRATION OVER ThL DUMMY VARIABLES P, Q, AND R. 

P , Q , R = CUMMY VARIABLES OF INTEGRATION 

PNT , GNT , RNT = CNE-DI MENSICNAL ARRAYS IN WHICH ARE STORED ThE 
VALUES OF THE ANSWERS TO THE INTEGRATIONS OVER THE DUMMY VARIABLES 
P, Q, AND R. 

R 3 = RATIO OF TUNNFL PERMEABILITY FACTTR TO BETA, R/BETA 

RP = A ONE-DIMENSIONAL ARRAY IN WHICH ARE STORED THE VALUES CF 
THE PERMEABILITY FACTOR OF THE TUNNEL 

SB = WING SEMI-SFAN TO TUNNEL SEMI-WIDTH RATIO, S /B 

SH = A CNE-DI MENSIONAL ARRAY IN WHICH ARE STORED THE VALUES CF 
THE RATIO OF WING SEMI-SPAN TC TUNNEL SEMI-HEIGHT, S/H. 

SMP , SMQ, SMR = CNE-DI MENS I CNAL ARRAYS IN WHICH THE INTERMEDIATE 
ANSWERS TO THE INTEGRATION PROCESSES ARE STORED. 

YH = A TWO-DIMENSIONAL ARRAY IN WHICH ARE STORED THE VALUES OF 
THE RATIO OF THE DISTANCE ALONG THE WING SEMI-SPAN TO THE TUNNEL 
SEMI-HEIGHT, Y/H. 

YS = NORMALIZED DISTANCE ALCNC- WING SPAN, Y/S 
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FORTRAN PROGRAM FOR CALCULATING UPwASH INTERFERENCE FACTORS 
AT THE CFNTtR OF WINGS CENTER-MOUNTED IN RECTANGULAR COMPLETELY 
CLOSED (ECUATION (3) IN MAIN BCDY OF REPORT) AND COMPLETELY OPEN 
(EQUATION (4) IN MAIN BODY CF REPORT) TEST SECTIONS. THESE 
EQUATIONS ARE INFINITE SUMMATIONS DERIVED BY THE METHOD OF IMAGES 
AS IN REFERENCE 2. 


IN ADDITION TC THE CCMMCN VARIABLE NAMES ALREADY CEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM. 

OL TC = CALCULATED LPWASH INTERFERENCE FACTOR FOR CLOSED TUNNEL 

DL TO = CALCULATED UPWASH INTERFERENCE FACTOR FOR OPEN TUNNEL 

I = INDEX ON VARIABLE SH 

J = INDEX ON VARIAELE RH 

K = INDEX ON SUMMATION. K IS THE SAME AS THE N OF EQUATIONS (3) 
AND (A). 

KCO = THE VALUE CF K AT WHICH THE SUMMATION IS TRUNCATED. 

ALL OTHER VARIABLE NAMES, EXCEPT AS NOTED, ARE USEO MERELY AS 
RQCKKEF P INC DEVICES. 


NAMELIST INPUT S ARE Eh, Sh 


PROGRAM*, T ( INPUT, CUT PUT ,TAPE5 = INFLT ,T AP E6 = CUT PUT , TAPE 1 , T AP E2 ) (B 1 ) 

DIMENSION BH(5) , SE< 2HNAMEL I S T / I NP T/BH , S B JR E AD< 5 , INPT) (8 2) 

PI=2.*ASIN(J.)*PI2=PI/2. $001 1=1 , 2$DClJ=i ,5$PBH=PI2*BH(J) (B3) 

A:. = PBH*SB(im:=i./(A.*Al*SB( I ) ) t T 2 = 1 . / ( 4. * SB ( I ) ) $TC1=-T2/SINH(A1) (B 4) 

TO 3=-T2/TANH( A. ) $A2=A3=0. *D02K = 1 ,25iARl=PEF*(2.*K+SB(I>) (8 5) 

AR? = PBH*(2.*K-SB ( I ) ) $ T4= 1 . / S I N M AR2 > - i . / S INH ( AR 1 ) $ A2 = A 2 + T4 (8 6) 

T5=< (-1. )**K)*(1./TANH(AR2)-I./TANH(AR1) )$A3=A3+T5 (8 7) 


T 4 AND T 5 RF PRESENT TFF SUMMATION TERMS IN EQUATIONS (3) AND (4), 
RESPECTIVELY, OF THE MAIN BCDY OF THE REPORT. WHEN BOTH T4 
AND T5 AP F LESS THAN O.OCOC1, THE SERIES ARE ASSUMED TO HAVE 
CONVERGED AND THE SUMMATION PPCCESS IS TERMINATED. 


IF ( ABS( TM .LT..CCO ) 1 . AND. AB S ( T 5 ) . L T . . 0 OOC'i ) G C T03 $ I F ( K . EQ. 25 ) GOTO 3 (8 8) 

2 CONTINUE (B 91 

3 KCO=K*TC2 = T2*A2$TUi = T2*A3tDLTC = U + TC] + TC 2 $D LTO=T H- TGi + T C 2 (B 10) 

WRITE (6,* ) SB ( I ) , Bh( J) , KCO , OL TC , D L T 0 (B 11) 

wR ITC ( 2,4) SB ( I ) , BH( J) , KC C , DL TC , D L T 0 (R 12) 

4 FORMAT (//2X*S/8=*FS.2,3X* ! B/H = *F6«2»3X*'SUMMATION INTEGER=*I 3»//3X*D (3 13) 

1BLTA, CENTER , CLOSED TLNNEL = *FR.4 , 5X*DELTA, CENTER, OPEN TUNN EL = *F 8 . 4 ) (B 14) 

1 CONTI NUESSTCPtCND (B 15) 

$ I NPT BH=o5,»75,l.,t.5,2.,SB=.?,.7,t (B 16) 
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FORTRAN PROGRAM FOR EVALUATING THE CLOSED-TUNNEL PORTION OF 
EQUATION ( A-2 ) OF APPENDIX A (THIS IS THE TERM CONTAINING TEE 
EXPONENTIAL FACTOR). THE RESULTS GIVE UPWASH INTERFERENCE FACTORS 
FOR WINGS CcNTEK-MOUNTFC BETWEEN INFINITE, VERTICAL, CLOSED 
BOUNDARIES. 


IN ACDITICN TC THE COM MTN VARIA8LE NAMES ALREADY CEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM. 

DLTC2 = THE UPWASH INTERFERENCE FACTOR CALCULATED FOR THE CLOSED- 
TUNNEL PORTION OF EQUATION (A-2), WHICH REPRESENTS CELTA-SUB-2. 

ALL OTHER VARIABLE NAMES, EXCEPT AS NOTEC, ARE USED MERELY AS 
BOCKKE E P ING DEVICES. 


NAMELIST INPUTS ARE EH 


PRUGR AMwT ( INPUT, OUTPUT ,TAPE5=INFLT ,T AP E6 =CUT PUT , T AP El , T AP E2 , PUNCH ) (B 17) 
CCMMCNSH (2 ) , I SH, YH ( 3 ,2 ) , IYF, BE, IBP (B 18) 

QI M C N S I ONFQ ( 3 ) , S MG (2 ) , ON T ( 1 ) , EH ( 5 > tNAMEL I ST / I NPT/ BH $RE AD ( 5 , INPT) (B 19) 

NQ = 2 $NFQ= 1 $JQ = 90$P 1 = 2. 3 A I 5 92 6 5 3 6 $ CNC = P I / 4 . $CNV=1. F-07tNMY=3 (B 20) 

THE FOLLOWING NEST OF CC-LCOPS INDEXES ON BH , SH, AND YH, 
CALCULATING THE VALUES OF SB, SH, YH, AND YS THEN USING THESE 
THROUGH THE SUBROUTINE SUBPROGRAMS INFINTQ, FUNCQ, AND MGAUSC TO 
FIND THE INTEGRAL AND, FRCM THAT, CLTC2. THE RESULTING DATA ARE 
PRINTED CUT AND PUNCHED ONTO C AT A- PROCE S S I NG CARDS FOP LATER 
COLLAT ICN IN ANOTHER PROGRAM WITH THE CATA REPRESENTING THE 
PROGRESSIVELY MORE OPEN TUNNEL. 


001 IRH=i, C $SH( _ ) =« 3#BH( I BH) $ SH (2 ) = • 7*8H( IBH) $DG1 I SH = 1 ,2 (B 21 ) 

SB = SH( I SH I / EH ( IBH) $BS = l./SBtWiR ITE ( 6 , 2 )BH( I EH ) »SB (B 22) 

WRITE < 2,2) BH< IBH) , Se (B 23) 

DO). I YH=1 ,NMYtYH( I Y H , I S H ) = S H ( l S H ) *F LO AT ( I Yh- 1 ) / FLO A T ( NMY- 1 ) (B 24) 

YS=YH(IYH,ISH)/SH(ISH)$CALLINFINTC(NC,SMG,FG,NFQ,JQ,CMC,CNV,QNT) (B 25) 

DLTC? = BS*QNT / (2 . *P I ) $W R I T E ( 6 , 3 ) Y S , DL TC 2 $ WR I TE ( 2 , 3 > Y S , DL TC 2 (B 26) 

2 FORMAT! 1 H ,///6X4riB/H=F7.3,lCX4HS/B=F7.3,//5X3HY/S,6X*0ELTA SUB C2 (B 27) 

1 *//) (B 28) 

3 FORMAT ( 1H , 2 XFt . 2 , 3 XF 10. 6 ) (B 29) 

PUNCH4,BH( IPH), SE»YS,CLTC2 (B 30) 

4 FORMA T ( 2 ( F 7. 3 ) » F 6. 2 » Flu • 6 ) (B31) 

1 CONT I NUEHSTOPtENC (B 32) 
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SUBROUTINE FONCQ IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST OVER THE VARIABLE C IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINE MGAUSC. ITS CALLING ARGUMENTS ARE Q, THE VARIABLE OF 
INTEGRATION WhCSE VALUE IS SUPPLIED BY MGAUSQ, AND FQ, THE NAPE OF 
AN ARRAY THAT WILL CONTAIN VALUES OF TEE INTEGRANDS COMPUTED 
DURING TEE INTEGRATION. 


SUBROUT INEFUNCQ ( C, FC ) JDIMENSICNFC 1 1 > , BH ( 5 ) (B 331 

CGMMONSHI2 ) ,ISH,YH( 3,2 )* IYH, BE, IBE (B 34) 

SHG=SH< ISH)*Q$YHC=YE( IYH, I SH ) * C $BHC = BH ( I 8H) *C (B 35) 

FQ = SINH(SHQ)*COSH<YH.G>/<EXP(f!HC)*S INH(BHG)) IRE TURN SEND IB 36) 


SUBROUTINE INFINTC IS A GENERALIZED SUBROUTINE SUBPROGRAM 
WRITTEN TO EVALUATE ^INFINITE* INTEGRALS BY MAKING USE OF SUB- 
ROUTINE MGAUSO (WHICH EMPLOYS THE GAUSS CUADRATURE METHOD IN INTEGRAL 
EVALUATION) IN SUCCESSIVE STEPS BEGINNING AT THE LOWER LIMIT OF 
INTEGRATION. THE INTEGRATION INTERVAL FOR EACH SUCCESSIVE INTEGRATION 
STEP IS SET EXTERNALLY BY TEE VARIABLE QINC IN THE CALLING ARGUMENTS 
OF INFINTC. THE INTEGRAL VALLE IS INITIALIZED AT ZERO ANC THE VALUE 
OF EACH SUCCESSIVE INTEGRATION STEP IS ADDED TO IT. THE INTEGER 
NFQ DETERMINES TEE NUMBER CF INTEGRANCS TC BE EVALUATED, WITH 
ANSWERS STORED IN THE ARRAY QINT. WHEN THE SUMMED VALUES OF THE 
INTEGRALS OVER SCME PARTICULAR STEP ARE LESS THAN CONV, THEY ARE 
ALL ASSUMEO TO HAVE CONVERGED AND THE INTEGRATION PROCESS IS 
TERMINATED. IF TEE NUMBER CF INTEGRATION STEPS IS LARGER THAN THE 
INTEGER JC, SET EXTFRNALL Y , TEE INTEGRATION PROCESS IS TERMINATED 
AND AN ERROR MESSAGE TEAT TEE INTEGRAL IS NON-CON VERG EN T WITH THE 
GIVEN UPPER LIMIT IS PRINTEC. THE CALLING PROGRAM MUST DIMENSION 
THE ARRAYS S U m Q , FCFQ , AND CINT AT THE MAXIMUM NUMBER TO BE USED 
IN THAT PROGRAM. 


SUBROUTINE I NF I N TC ( NC , SUM C , FC F C , NFC , J C, C I NC , CON V , C I NT ) (8 37) 

DIMENSION SUMQ( i ) , FCFQI 1 I ,QINT( 1 ) $0011 = 1, NFC (B 38) 

1 TINT! 1)='.' o$AC=C.I0021Q=l,JG$IF(IC-l)4»4,3 (B 39) 

3 AC = BO (B 40 ) 

4 RQ=FLGAT( I C ) *CI NCSC AL LMG AUSQ < A C , EC , NC » SU MG , FOFQ , NFC ) (B 41) 

<1CUNV=R. $D05I=.» , NF0ICC0NV = CCCNV + AES(SUMQ( I) ) (B 42) 

5 CINT! I)=GINT< I ) + SUM C I I) $ I F (QCC NV . LT . CONV ) GOTO 6$ I F I I Q . EQ . JQ )GCT07 (B 43) 

G0TO2 (B 44) 

7 PR I NT 8 $GC T 06 (8 45) 

R FORMAT! 1H ,2XMT NTFGRAL NQNCCN V E RC EN T WITH GIVEN UPPER LIMIT ON C*> (B 46) 

2 CONTINUE (B 47) 

6 Kt TURNttNU (6 48) 
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THE: SUBROUT I Nfc MGAUSQ EMPLOYS THE GAUSS QUADRATURE METHCD TO 
EVALUATE * NUMBER * NUMBER CF INTEGRANDS F(X)DX BETWEEN THE LIMITS A 
AND B. ANSWERS ARE STORED IN THE ARRAY SUM. FUNCQ IS THE NAME OF A 
SUBROUTINE SUBPROGRAM USED TO EVALUATE ThE INTEGRANDS. THE NUMBER 
OF POINTS USEC WITHIN ThE INTEGRATION LIMITS IS TEN TIMES THE 
INTEGER N. THIS SUBROUTINE SUBPROGRAM IS A PRELIMINARY VERSION OF 
THE SUBROUTINE MGAUSS NOW ON THE LIBRARY TAPE OF THE CDC COMPUTER 
SYSTEM AT LANGLEY RESEARCH CENTER. A NCR E COMPLETF CISCUSSION OF 
IT THAN GIVEN HERE IS GIVEN IS SECTION D 1.1, VOLUME I, OF THE 
LANGLEY RESEARCH CENTER COMPUTER PROGRAMMING MANUAL. 


SUBROUTINE MG AU SC t A , B , N , SUM , FCFX , NUMB ER ) (B 49 » 

DIMENSION U(5), RI5), SUM! 1 ) , FGF X ( 1 ) SDO ILL = 1 , NUMBER (B 50) 

1 SUM! LL ) =Co (B 51) 

IF< A. EQ.R) RETURN $U( 1 (=.4 25 56 26 30 51,51 64 $U( 2 )=. 283 302 3 02 9 8 5376 (B 52) 

U<2)=. 1 6f 29521555 94 66$U<4)=.G6746 83166 55 50 8 $U(5)=. 9 130 46735741414 (B 53) 

R (1) = . 1477 6 2 li 23 5 73 76 SRI 2 ) =. 1 2 4 6 3 33 5 56 5495 8 $R ( 3 ) =. 1C 9 543 1 812 5799 1 (B 54) 

R(4)=. 9 747256 74 57529CFR(5)=. 0333356721 54 344$FINE=N (B 55) 

OELTA = F INE/I B-A ) $D0 3K =1 , N $X I = K- 1 $F I NE= A+ X I /DEL T A$ 002 1 1 = 1 ,5 ( B 56 ) 

UU=U(II)/DELTA+FINE$CALLFUNCC(UU,FCFX) $D02JCYB0Y=1, NUMBER ( B 57) 

2 SUM ( JO YBOY ) = R ( I I ) 4FQF X( JOYBOY ) + S L M I JO YB CY ) $D03 J J = 1 , 5 (B 58) 

UU=(i.9-U< JJ) )/0ELTA+FINE$CALLFUNCQ(UU,F0FX I tD03NN= 1 , NUMBER (B 59) 

3 SUM! NN ) =R ( JJ ) *FOFX( NN )+SLM(NN) SDC7I JK=1, NUMBER <B 69) 

7 SUM ( IJK) =SUM( IJK )/CELT A $ RETU F N $ EN C (B 61) 

$ I NPT BH=.5,.75,3 .,1.5,2., $ (B 62) 
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FORTRAN PRCGRAM FOR EVALUATING THE C LC SED- TUNNE L PORTION OF 
EQUATION (A-4) OF APPENDIX A (THIS IS THE TERM CONTAINING THE 
SINGLE INTFGRAL). THE RESULTS GIVE UPW4SH INTERFERENCE FACTORS 
FOR WINGS CENTER- MOUNTED BETWEEN INFINITE, HORIZONTAL, CLOSED 
BOUNDARIES* 


IN ADDITION TO THE COMMCN VARIABLE NAMES ALREADY CEFINEO, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TO THIS PARTICULAR PROGRAM. 

ARGYS = THE VARIABLE MULTIPLIED CNTO THF VARIABLE OF INTEGRATION 
TO MAKF UP THE ARGUMENT OF THE SINE TERN CF THE INTEGRAND. 

DL = THd VALUE CF THE I NT E 0 P A l RESULTING FROM USE CF A POSITIVE 
SIGN IN THE VARIABLE ARC-YS. 

02 = THE VALUE CF THE INTEGRAL RESULTING FROM USE OF A NEGATIVE 
SIGN IN THF VARIABLE ARGYS. 

QLTC3 = THE UPNASH INTERFERENCE FACTOR CALCULATED FOR THE CLOSED- 
TUNNEL PORTION OF EQUATION I A-4 ) , WHICH REPRESENTS DELTA-SUB-3. 

ALL OTHER VARIABLE NAMFS, EXCEPT AS NOTED, ARE USED MERELY AS 
BOOKKEEPING DEVICES. 


NAVEL 1ST INPLTS ARE BH 


PROGR AMWT ( INPUT , OUT PUT , TAPE5 = lNPLT,TAPEc> = OUTPUT, TAPE I , TAPE? , PUNCH) (B 63) 

01 MFNSIONFQO. ) , S MG < 1 ) , QN T ( 1 ) , S H ( 2 ) , YH U , 2 ) , BH ( 5 ) ( 8 64 ) 

COM MON ARGYS* NAM EL I S 1 / IMP T/BH tRE AC ( 5 , I N PT ) (B 65) 

NQ = 5 1 NF Q = \ TJC=6 ‘ *F I =3. i 41 592 6 5 36 tCNV=l .E-07 *NMY = 3 < 3 66) 

IN' ADDITION TC THE OPERATIONS PERFORMED BY THE NESTED DO-LOOPS 
OF THE D R E Ct DING PROGRAM, THOSE CF THIS PROGRAM PERFORM TwC SEPARATE 
INTEGRATIONS. IN THE FIRST, THE VARIABLE MULTIPLYING THE INTEGRATION 
VARIABLE IN THE ARGUMENT DF TEE SINE TERM OF THE INTECRAND HAS A 
POSITIVE SIGN. IN THE SECOND, IT HAS A NEGATIVE SIGN. FOR EACH 
INTEGRATION, CMC -- THE VARIABLE SETTING THE INTEGRATION INTERVAL -- 
IS A FUNCTION CF ARGYS. THIS PROCEDURE RESTRICTS THE INTEGRATION 
LIMITS CN ANY PARTICULAR INTEGRATION INTERVAL SUCH THAT THE 
ARGUMENT CF THE SINE TERM IN TEE INTEGRAND IS SOME MULTIPLE CF 
PI AT THE LIMITS. 


DCiI8H=l , 6BSHM )=.^*fiMIPM*SM2) = .7*eH I BH ) *nOiISH=l, 2 ( B 67) 

5>B=SH(ISH)/eH<I3H)$BS = l./S6SnRITE(6,2)3H<IPHI,S8 ( B 68 ) 

WRI TE (2,2 ) BHIIBH ) , SB (B 69) 

2 FORMA T U H ,///b<4HB/H=F7.3,i:X4HS/e = F7o3,//SXbHY/S,8X2HC] ,i2X2HD2, (B 70) 
13X*DELTA S’JH C?*//) (B 71) 

DOE I YH = 1 , N my $YH ( IYH,ISH) = SH( I SH) *FLOAT ( I YH-1 ) /FLOAT ( NMY-U {B 72) 

AKGYS=YH( IYH,ISH)+SH( ISH) * Y$ = YF ( IYH, ISM /SH( ISH) (B 73) 

QNC=PT/APGYS$CALLI NF INTQ( NC, SRC, FC , NF C , J Q , C NC , CN V , ; NT ) 40 J =QNT <8 74) 
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IF IYH = N MY t THE VARIABLE ARCYS = J ANC THE VALUE CF THE INTEGRAL 
IS ZFRP. 

IF( I YH„ EQ.NMY )GCTC6 (B 75) 

ARGYS = YH< I YH, ISH )-SF< l SH ) $GNC = P I / A f S ( A RG Y S J $YS=YH ( I Y I- , ISP )/ SHI ISH) ( S 75) 

C A L L I N F I NT G I N C * S MG , FC * NFQ , JQ , C NC ,CNV , QNT ) 1D2=QNT$G0T07 (B 77) 

6 02 = ). (B 78) 

7 r>LTC3=8S*( C1-D2 ) /<4.*Pl) $WRI TE (6 , 3 ) YS , Q1 . C2 . DLTC 3 (B 79) 

WRITE(2,3) YS,0i t 02,CLTC3 IB 80) 

3 FORMAT ( i H , 2 X Fb . 2 ♦ 3 ( 4X FI 0 . 6 ) ) <B 81) 

PUNCH4,BH< IBM) , SB, YS,DLTC3 (B 82) 

4 FURMAT (2 ! F 7 .3I, F 6 » 2 » F 1 . 6 ) (B 83) 

1 CONTINUE 1ST0PTENC (B 84) 

SUBROUTINE FLNCC IS CNE WRITTEN TO EVALUATE THE INTEGRANC OF 
INTEREST OVfcR TEE VARIABLE G IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINE MGAUSC. REFER TC CISCUSSICN IN THE PRECEDING PROGRAM. 

SUBROUTINE FUNCCi ( G . FQ ) $ CCMMC N AFGY S t A RG = ARGYS*Q (B 85) 

F G = S IN I ARG)/ IEXP (0 ) *COSH ( Q) ) $RETLRN$END (B 86) 

$INPT 8H=.5,.75,J .,1.5,2. $ (B 112) 

THIS PROGRAM MAKES USE CF THE SUBROUTINE SUBPROGRAMS INF I NT Q AND 
MOAUSC AS DEFINLC AND CISCUSSEC IN THE PRECEDING PROGRAM. 
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FORTRAN PROGRAM FOR EVALUATING EQUATICN (A-7) — WHICH REP- 
RESENTS THE UPWASri INTERFERENCE FACTORS OUE TO THE EFFECT OF 
HORIZONTAL CLCSEC BOUNDARIES CN THE INTERFERENCE VELOCITY POTENTIAL 
INSIDE THE TEST SECTION DUE TC VERTICAL CLOSED BOUNDARIES — OF 
APPENDIX A ANC EQUATION (A-8) — WHICH REPRESENTS THE EFFECT OF 
HORIZONTAL CLCSED BOUNDARIES CN THE INTERFERENCE VELOCITY POTENTIAL 
OUTSIDE THE TEST SECTION DUE TC THE HORIZONTAL CLOSED BOUNDARIES. 
THIS PROGRAM CIFFERS FR CM THE PRECEDING CNES IN THAT COUBLE IN- 
FINITE INTEGRALS IN TWO VARIABLES ARE CALCULATED RATHER THAN THE 
SINGLE INFINITE INTEGRAL IN ONE VARIABLE CALCULATED IN THE PRE- 
CEDING PROGRAMS. 


IN AUDIT I CN TO THE COMMON VARIABLE NAMES ALREADY DEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM, 

Ai = A ONE-DIMENSIONAL ARRAY IN WHICH APE STORFC THE VALUES OF THE 
INTEGRALS RESULTING FROM USE CF A POSITIVE SIGN IN THE VARIABLE *BY*. 

42 = A ONE-DIMENSIONAL ARRAY IN WHICH ARE STORED THE VALUES CF THE 
INTEGRALS RESILTING FROM USE CF A NEGATIVE SIGN IN THE VARIABLE *BY*. 

BY = A VARIABLE WHICH SERVES THE SAME FUNCTION AS THE VARIABLE 
ARGYS IN THE PRECEDING PROGRAM. 

DLTC4 = THE UPWASH INTERFERENCE FACTOR CALCULATED FCR EQUATION (A-7). 

DLTC5 = THE UPWASH INTERFERENCE FACTOR CALCULATED FOR EQUATICN (A-3). 

ALL OTHER VAR IABLfc NAMES, EXCEPT AS NOTED, ARE USED MERELY AS 
BOOKKEEPING DEVICES. 


NAMEL 1ST INPLTS ARE BH 


PROGRAMWT ( INPUT, CUT PL T»TAPE5=INPLT»TAPE6 = CUTPUT»TAPE1»TAPE2» PUNCH) ( R 113) 

CUMMONSH12 ) , I S, YH< J ,2) , IY.BH, IE,C,BY, PI (B 114) 

DIMENSION A 1 ( 2 ) , A2 ( 2 ) (B 115) 

DI MENS I ONFQ ( 2 ), S MQ( 2 ), QNT ( 2 ) »BH(5)$NAMELIST/INPT/BH$READ(5» INPT ) (8 116) 

NQ=3$NFQ = 2*J« = 5'J$PI = 2.1^i 592£ 536*CNV = i.E-‘J7*NMY = 3 (B 117) 

DOlI8=l,5$SH(l)=.3*EF< IB)$SH(2 )=.7*0H( I B ) $DC 1 1 S= 1 , 2 (B 118) 

SB=SH( IS) /0H< IB ) *BS=1. /S6$WRI TE ( £ ,2) BH ( I B) , SB (B 119) 

WRITE (2 ,2 ) BH( IB) ,SB ( B 120) 

2 FORMAT ( 1H , / / /6 X4HB /h = F7 . 3 , 10 X4H S /B = F 7. 3 , / / 5 X3HY /S , 7 X*DEL T A SUB C4 (B 121) 

1*,2X*DELTA SUB C5*//) (B 122) 

D01I Y = l,NMYJYM( I Y, I S» = SH( I S> *FLOAT ( I Y-H /FLOAT (NMY-1 ) ( B 123 ) 

8Y=HH ( IB) +YF( IY , IS) $GNC= P I (B 124) 

YS = YH( IY,IS)/SH(IS) tCALL I NF I NTQ ( N Q , SMQ , F Q , NFQ , J Q , QNC ,CNV,QNT ) (B 125) 

Aid) =QNT (l)tAj (2) = CNT(2)tBY = BH( IE) — YHI IY, IS) (B 126) 

CALLINFINTQ< NQ, SMQ , FC ,NFQ, JQ ,QNC ,CNV,CNT ) $A2 ( 1 )=QNT( 1 ) (B 127) 

A 2 (. 2 )=QNT( 2 ) inLTC4=-eS*<Al(i)4A2(l) )/(2.*PI*PI ) (B 123) 

0LTC5=BS*< Al ( 2) + A2 ( 2) > /(2.*PI*P I )$WRITE( 6 ,3 ) YS , DLTC4 , DL TC 5 (B 129) 

WRITE(2,2)YS,DLTC4,CLTC5 (B 130) 

3 F0RMATI1H , 2 XF6 „ 2 , 2 ( 4 X FI . 6 ) I (8 131) 

PUNCH4,BH( IB),SB,YS,DLTC4,DLTC5 (B 132) 

4 FOPMAT(2(F7.3),Ffc.2 , 2 ( F 3. J . 6 ) ) (B 133) 

1 CONTI NUFtSTCPTENC (B 134) 
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IN ACCITICN TC PERFORM IMG THE FUNCTIONS PERFORMED BY SUBROUTINES 
OF THE SAME NAME IN PRECEDING PROGRAMS, THE SUBROUTINE FUNCO IN THIS 
PROGRAM MAKES USE OF THE SUBROUTINES INFINTR AND MG A U SR IN ORDER 
TO PERFORM THE SECOND INTEGRATION OVER THE VARIABLE R REQUIRED BY 
EQUATIONS IA-7) AND (A-8). 

SUBROUTINE FUNC Q ( Q , F Q ) $D I ME N S I ONFQ ( 2 ) , SMR I 2 ) , FR ( 2 ) , RNT ( 2 ) , BH ( 5 ) (B 135) 
COMMON SHI 2 ) , I S , YH ( 3 ,2 ) , I Y , BH , I B , S , BY , P I $S = C $NR = 2 $NFR= 2 t J P= 50 ( B 136) 

IN THIS FROGRAM, THE INTEGRATION INTERVAL RNC IS A FUNCTION OF 
THE VARIABLE *BY * ANALOGOUSLY TC AND FOR THE SAME REASONS THAT THE 
INTEGRATION INTERVAL QNC IS A FUNCTION OF THE VARIABLE ARGYS IN THE 
PRECEDING PROGRAM FOR DLTC3. 

RNC=P I / ABS ( BY ) tCNV=! . F-07JC ALL INFINTR (NR ,SMR,FR,NFR,JR,RNC,CNV,FQ) IB 137) 
RETURNtEND <B 138) 

SUBROUTINE FLNCR IS CNF. WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST 0 VCR THE VARIABLES C AND R IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINES INFINTQ, INFINTR, MGAUSG, AND MGAUSR, 


SUBROUTINE F UNC R I A , F R ) $ C I ME N S IC N F P 1 2 ) , BH I 5 ) IB 137) 
CO TMONSH I 2 >,IS,YF<3,2),IY,BH,IE,C,BY,PI$SHC = SH(IS>*Q$BFG = BH(IB)*Q IB 140 ) 
ARG1=C0S(Q)*SINH(SHC)/ I E X P I BHC )* I C *Q + R *R )*COSHIR) )$BYR = BY*R IB 141) 
ARG2=(Q*S INHI BHQ ) *C C S I BYR ) +R*CC SH I BHC ) *SIN(BYR))/SINH(BHC) IB 142) 
ARG3=C*COS(fiYR>-K*SIN(BYR)$FR<l>=ARGl*ARG2$FR(2)=AR61*ARG3 IB 143) 
RETURNtEND IB 144) 
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THE SUBROUTINE SUBPROGRAMS INFINTR AND MGAUSR IN THIS PROGRAM 
ARE NECESSARY TC ALLOW INTEGRATIONS OVER TWO VARIABLES RATHER THAN 
THE ONE VARIABLE OF PRECECING PROGRAMS. EXCEPTING FOR THE SUBROUTINE 
NAMES ANC WHATEVER VARIABLE NAME-CHANGES ARE NECESSARY TO ALLOW THE 
CALCULATIONS TC PROCEED WI1FOLT CONFUSION, THEY ARE IDENTICAL WITH 
THE SUBROUTINES INFINTG AND MGAUSC AS DISCUSSED IN A PRECEDING 

PROGRAM. 


SUBROUTINE I NF I NTR ( NS , SUM G , F OF R , NF Q , JQ , C I NC , C ON V , Q I NT ) 

DIMENSION SLMQIl ),FCFR(i),CINT(l) $DC1I = 1,NFQ 

1 >J I NT( I ) =>_ » tAO='. » $00 210=1 , JQi IF( I C — I. >4,4,3 

3 A C = B Q 

4 3Q = F LO A T ( I Q) *QI NCtCALL MGAUSR (AC, B C ,NC , S U MC » FCFR, NFQ) 

QCCNV = ).$D05 I = i , NFiJ$CCONV = QCCNV + ABS{ SLMQ ( I ) ) 

5 QINT( I) = 0 I NT I I ) + SUMC( I ) i IF ( QCC N V . L T . CC NV 1G0TC6SI F ( IQ.EQ. JQ)GCT07 

GOT 02 

7 WRI IF (6,8) 5WRITL (? ,8) SG0TC6 

8 FORMAT ( 1H ,2X*INTEGRAL NONCON V ERG E NT KITH GIVEN UPPER LIMIT ON R*I 

2 CONTINUE 

6 RETURNSENO 


SUBROUTINE MGAUSR ( A , E , N , S UM , F C E X , NUM B ER ) 

DIMENSION i) <5 ),R <5 ) ,SUM( i ) ,FOFX< ] ) t DO ILL = I , NLMBE R 
1 SU V ( LL )=u 0 

IF ( A.EQ.B (RETURN $U( 1 )=„425 56 2 6 30 5}9184tU< 2)=. 2 83 3>" 2 3C2 ° 3 53 76 
U( 3) =„ 16 r 2 05 21 5 5 5 :,4E6I>U< 4 > = .' 67468316655 50 8 $U ( 5 ) = . 11 30 467 357 41 A 1 A 
R< i )=. 1477 421) 235737e*k! 2 ) =. 1 2 A 6 2 2 3 5 c . 6 5 454 £ $ P ( 3 I =. 1 0 <5 543 1 81 2 5799 1 
<<< A) =. 1747 2 5 67^5 752 90*R< 5 »■=. t-3 3 3 2 5 67 2 1 54 3M $FINE=N 
DELTA=F INF/( 0-A ) lD0 3K = j ,NtXI =K-1$F I NE = A+ X I /DEL TA$DC2 I 1 = 1 ,5 
UU = U( I I I/UELTA+F INE tCALLFU NCR (UU,FOFX> $002JQYR0Y = J. , NUMBER 
? SJM( JCYBOY>=R(I I )*FCFX< JOYBOY 1+ S LM ( JC Y0 C Y ) $ DO 3 J J = 1 , 5 

UU = (i® : -U ( JJ ) ) / 0 ELT A + FINF$CALLFUNCP (UU,FOFX ) $D03NN= 1 , NUMBER 
3 SUM(NN) =P ( JJ )*FUFX< NNI+SUMCNN) I DC 7 I JK = L, NUMBER 
7 SUM( I JK) =SUM( IJKI/DfcLTA SRETURNtENC 

$ I NP T *3 H — « *5 j#7? f I# * F * 2 © » $ 


(6 145) 
(3 146) 
( 0 147) 
(B 148) 
( B 149) 
(B 150) 
(B 151) 
(B 152) 
(B 153) 
(B 154) 
(B 155) 
(B 156) 


( B 182) 
(B 183) 
( (3 184) 
(B 185) 
(B 186) 
(B 187) 
( B ] 88) 
(B 189) 
(B 190) 
(B 191) 
(B 192) 
(B 193) 
(B 194) 

(8 195) 


THIS PROGRAM M^KFS USE OF THE SUBROUTINE SUBPROGRAMS INF I NT Q ANO 
MG AU SO AS HEFINEC ANO CISCLSSEC IN A PRECOOING PROGRAM* 
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FORTRAN PROGRAM FOR COLLATING DATA- THE PUNCHED-CARO OLTPUT 
OF THE THREE PRECEDING PROGRAMS — WHICH REPRESENTS UFWASH INTER- 
FERENCE FACTORS FOR WINGS C EN T ER-MOUN TF D IN CLOSED RECTANGULAR WIND 
TUNNELS AS DEFINED IN THOSE PROGRAMS ANC APPENDIX A — IS THE INPUT 
TO THIS PROGRAM. THE CATA AS READ INTO THE PROGRAM APE PRINTED- 
OUT, THEN COLLATED AND RE-PRINTED AS A CHECK. THE COLLATED DATA 
ARC ALSO CUT PUT CN PUNCHED C A T A-PRCCE S S I NG CARDS FOR LATER COLLATION 
WITH THE DATA REPRESENTING UPWASH INTERFERENCE FACTORS FOR WINGS 
CENTER-MOUNTED IN P ROGRF S S I V E L Y MORE OPEN RECTANGULAR PERFORATED 
WIND TUNNELS. 


IN ADDITION TC THE COMMON VARIABLE NAMES ALREADY DEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM. 

02 = DLTC2 OF A PRECEDING PROGRAM. 

DR = DLTC? OF A PRECEDING PROGRAM. 

DA = DLTC4 OF A PRECEDING PROGRAM. 

05 = DL TC f OF A PRECEDING FRCGRAM. 

ALL OTHFp VARIAbLF NAMES, EXCEPT AS NCTED , ARE USED MERELY AS 
BOCKKFFPING DEVICES. 


PROGR AMWT ( I NPUT, OUTPUT ,T APE 5= INPLT , TAP E6 = CL TPUT , TAPE 1 , TAPE2 , PUNCH) 
01 MFNSIUNRH (30,3 ) , S E < 3C , 3 ) , YS < 31 ,3 ) ,02 < 30 » , C3 < 30 » , 04 < 30 > , 05 < 30 ) 
READ! 5, 1 ) ( PHI 1,1 ) ,S8( I ,1 ) ,YS( I ,1) »C2( I ) , I = i,3'J) 

READi 5,1 ) < BH( I ,2 ) ,SE < I ,2 ) ,YS< I ,2 ) ,03 < I) , 1 = 1 , 3'') 

) FnRMAT(2(F7«3),Ffc.2,F5.0.6) 

RlAD< 5,2 ) < BH< I ,3 ) ,SE ( I ,3 ) ,YS( 1 ,3 ) ,C4 < I ) ,D5 < I ) , 1=1, 30 ) 

2 FORMAT (2 (P'7. 3), Ft. 2, 2 (FID. 6 I ) 

WK I T E < 6 , :• ) ( 3H < I , 1 I , SE < 1 , 3) , Y S ( I , H , 02 < I) , I = 1 , 30 ) 

WR ITE (6,f ) ( I3F< I , 1 1 , SR < I, 2 ) , YS< I , 2) ,03 < II , 1 = 1 , 30) 

WPI Tl(6,a) < BH (I , 3 ) , SB ( I ,3 ) ,YS< I ,3 ) ,D4( I I ,C5 (I I , 1=1 , 30) 

WR ITE (6,6 ) < SB< I , i ) ,BH< I , 1) , YS< I , 1) ,D2< I ) ,D3< I ) ,04(1) ,05 < I ) , I =1 ,3 3) 
WRITE<2, F)(SP<I,l),EH(I,i),YS(I,l), 02(1), C2(I)»D4(I), 05(1), 1=1, 30 

3 FORM AT ( 1 F 1 , //2X*S/0*4X*B/H*3X*Y/S*2X*DELTA SLB C2*//(2(F7#3) , F6. 2 , 
IF' 1 2. t) ) 

A FORMAT! 1H , / /2 X * S /R *4 X*B /H <<3 X* Y / S*2X*DEL T A SUB CA*4X*DELTA SUB C5* 
1//(2(F7.2),F6.2,F12.6,AXF)2.6)I 

5 FORMATtiH ,//2X*S/R*4X*B/H*3X*Y/S#2X*DE0TA SLB C3*/ / ( 2 < F 7 . 3 ) , F6. 2 , 
i FI 2.6 I ) 

6 FORMAT ( !H , / / ?X * S /B * A X *B /H* 3 X* V / 5* A X*DEL T A SUB C2*4X*DEUTA SUB C3* 
1 AX*DE U T A SUB C4*4 X*CE LT A SUB C 5 * / / ( 2 ( F 7 . 3 ) , Ffe . 2 , A ( AXF1 2 . 6 ) ) ) 

I A=-2 too! " J = 1 ,2* I A= I A+3$MI=I A- 6 $DC1CK = 1 , 5 $M= M+6 t N = M + 2 
WRITE(fc>,7)(S3(I»l)*EH(I*5.)»YS(I*3)»D2(l)*D2(I),04(I)»D?(I)»I=M,N) 
WRI TE ( 2,7) ( SB( I , 1) , BH < I , l ) ,YS( I ,JI ,02 ( I ) ,D3< I ) ,04 < I ) ,D5 ( I) , I =M,N) 

7 FORMAT (2 ( F7.3), F E • 2 , 4 ( 2XF12.6 ) ) 

PUNCFI8 , ( SB ( I , 1 ) , BH < I , 1 ) , YS < I , i I , C 2 < I ) , 03 ( I I , 04 ( l> , D5 < I ) , I = M , N > 

8 FORM AT (2 (F7. 3), Ft. 2, A (FI 2. 6) > 

1C CONTINUE$STOP*ENC 


( B 

196) 

( B 

197) 

( B 

198) 

( B 

199) 

( B 

200) 

( B 

201 ) 

(B 

202) 

( B 

203) 

( B 

204) 

( B 

208) 

( R 

206) 

( B 

207) 

( B 

208) 

( B 

209) 

( B 

210 ) 

( B 

211) 

( B 

212) 

(8 

213) 

(B 

214) 

( B 

215) 

(B 

216) 

( B 

217) 

( B 

218 ) 

( B 

219) 

( B 

220) 

( B 

221) 

( B 

222) 
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F n K THAN PROGRAM FCP EVALUATING THE PROGRESS IV ELY-MCRE-OPEN- 
TUNNEL PERT l UN DF ECU AT ION <A-2) CF APPENDIX A (THIS IS THE TERM 
CONTAINING THE HYPERBOLIC SECANT IN THE INTEGRAND),, THE RESULTS 
GIVE UPWASH INTERFERENCE FACTORS FOR WINGS CE NTER -MOUNTED BETWEEN 
INFINITE, VERTICAL, PERFORATED BOUNDARIES. 


IN AUDITION TO THE CCMNCN VARIABLE NAMES ALREADY CEFINED, 

THE FOLLCWING NAMES ARE OF IMPORTANCE TO THIS PARTICULAR PROGRAM. 

DLT 02 = THE IJPWASF INTERFERENCE FACTOR CALCULATED FOR THE PROGRESS- 
I VELY-MCPC-JPEN-PLRFORATEr-TUNNEL PORTION OF EQUATION { A- 2 ) . 

ALL CTHtP VARIABLE NAMES, EXCEFT AS NOTEC, ARE JSEO MERELY AS 
BOOKKt E P I \G DEVICES. 


NAMELIST INPUTS ARE FF, BT , FP 


PRf GRAMWTI INPUT, CUT FO T ,T A P E 5 = INPUT ,TAPE6=CLTPUT , TAPE1 , TA P F 2 , PUNCH) ( R 22 3) 


COMM )N SHI 2 ) , I S, YF( 3 ,2) ,1 v ,BH, IB , PI ,RP IB 2 24) 

H w cN S I ON F 0 ( '. ) , S W C ( ,. ) ,QNT (1) » BH( S ) , PT ( 5 ) ,P P ( 4 ) ( R 22 5) 

NAMbLIST/INPT/3i,HT,RP$READ(5, INPT ) (P 226) 

l\C = 4tNP£ = 1 S J C = T ) . *p 1=3. 14 169^6 S3£ JC NV= 1 • F -L 7 $NM Y= 3 (B 22 7) 

DO ’ l S = \ f J 1> J) r , 10 = *., 6 $SH ( 1) = 0 ?*0H( IB) $SH(2 )=.7*HH( IB ) $DO;. IP = 1,4 ( R 228 ) 

i;p'IT=l,J $RP = RP( IP )/BT { IT ) (B 220) 


IN THIS FFOGRAM, THE INTEGRATION INTERVAL QNC IS m A DE A FUNCTION 
OF TOP VARIABLE RE AND THE INDEXING INTEGER IP IN ORDER TO GIVE THF 
SAME NUMBER OF INTEGRATION PCINTS PER UNIT LENGTH OF THE INTEGRATION 
VARIABLE G IN THE ARGUMENTS CF THE TRANSCENDENTAL FUNCTIONS IN THE 
INTEGRAND. THIS ALSO PRODUCES APPROXIMATELY THE SAME AMOUNT CF 
COMPUTER PROCESSING TIME FCR TEE EVALUATION OF EACH TEST CASE. 


JNC=?R/Fl OAT ( IP ) *>S0 = SH( I S) /BH( IB ) *BS=1 . /SB (023^) 

WF IT,: ) Pm ( IH) t SB tRPfWRITE ( E ,2 )PH( IB ) , SB,RP ( B 2 3/) 

2 FOr-viATHh , ///F.X^M 5 /h=F7.3 , 10 X4HS / 3 = F7 . 3 , " sX*4/BETA= *F7. 3 , // 5X3FY/ {R 232) 
1 S,7X*CFLTA SUP 0/*//) (B 233) 

on: I Y=i ,N"Y$YH( I Y, I S) =SH ( I S) *FLCAT ( I Y-C ) /FLOAT (NMY-l ) ( P 234) 

YS = YH( I Y, I S ) /S U IS ) iCALL INF INTO! N C , SMC ,FQ,NFC, JQ ,QNC ,CNV ,QNT ) ( B 233) 

0LTO2 = BS*UNT/(2.*PI ) $WR I TF (6 ,3 ) Y S , DL TC2 $WR I TE ( 2, 3 ) Y S , IU TC2 ( 8 236) 

PUNCH4, SR, Rh( IB) ,KB ,YS,PLTC2 (B 237 ) 

4 FCPMATI 3(F7.?) ,F6. C ' ,Fi;.6) (B 238) 

A FORMAT ( 4 * F t> c ? , 6 X F 1 < ' # 6 ) (R 239) 

? CONTI NUE ^STCPiENC <8 240) 
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SUB F; CUT I NG EUNCC IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTERtST OVER THE VARIABLE 5 IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINE MGAUSCo REFER TC DISCUSSION IN A PRECEDING PROGRAM. 


SUBROUTINE- FUNCG(C.FO) $CIMENS IONFQ(i) , BH( 5 ) (R 241) 

CONMONSHI? ) iIStYHT »2) » I Y » RH , 16, FI , R B (B 24 2) 

S0 = SH(IS)*-nv3=YMIY»IS)*QSB0=eH(IB)*GtBR*l./RB (B 243) 

A1=SI NM( SG)*CrSM(YO)/I SINHIBC )*CCSH(HC ) » (B 244) 

A2= SORT (RE #BR+T ANh( BC ) *T ANH ( R G ) + I . ) ( 8 245) 

FQ = A1*(HH*TANH BCI-A2 ) 4RETURN$ENC IB 246) 

iINPT RH = . 5 ,. 75 ,1 . . 5,2 . ,BT = 3 » ,.8 , .6 , .45 ,.3, RP=.I , .45, 2. , 7. 5, $ ( 8 272) 


THIS PROGRAM MAKES USE OF THE SLBRCLTINE SUR D F GORANS INFINTC AND 
v GAUSu AS Of.FINEC AND CISCLSSEC IN A PRECEDING PROGRAM. 
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FORTRAN PROGRAM FOR EVALUATING THE PROGRESS I V ELY- MGR F -OPEN- 
TUNNEL PCRTIGN OF EQUATION (A-4) CF APPENDIX A (THIS IS THE DOUBLE- 
INTEGRAL PORTION OF THE EQUATION). THE RESULTS GIVE UPWASH INTER- 
FERENCE FACTORS FOR WINGS CENTER — MG UN TED BETWEEN INFINITE, HORIZONTAL, 
PERFORATED BOUNDARIES. 


IN ADO IT I CN TC THE CCMMCN VARIABLE NAMES ALREADY DEFINED, 

THE FOLLOWING NAMES ARE OF IMFCRTANCE TC THIS PARTICULAR PROGRAM. 

AYS * A VARIABLE *HICH SERVES THE SAME FUNCTION AS THE VARIABLE 
ARGYS IN A PRECEDING PRCGRAM. 

Qi = THE VALUE CF THE INTEGRAL RESULTING FROM USE OF A POSITIVE 
SIGN IN THF VARIABLE AYS. 

D2 = THE VALUE CF THE INTEGRAL RESULTING FROM USE OF A NEGATIVE 
SIGN IN THE VARIABLE AYS. 

DLTQ3 = THF UPWASH INTERFERENCE FACTOR CALCULATED FOR THE PROGRESS- 
IVE LY-MCRE-OP EN-PERFOR A TEC-TUNNEL PCRTICN OF EQUATION (A-4). 

ALL CTHEF VARIABLE NAMES, EXCEPT AS NOTEC, ARE USED MERELY AS 
BOOKKEFP INC DEVICES. 


NAMELIST INPLTS ARE p F , FT, PP 


P K OG R A M W T ( INPUT, CUT PUT ,T APE5= INPLT , T A P E6 =CUT PUT , T AP E 1 , T A P E 2 , PUNCH ) { B 273) 

COMMONSH { 2 ),IS,YF(3,2),IY,BH,IG,C,PI,R8,AYS,IP (R 274) 

DIMENSIONED*;, ) ,SMQ(1) , QNT ( 1 ) , BH ( 5 ) , BT ( 5 ) ,RP(4) (B 275) 

NAMfcL I ST/INPT/BH ,RT,RP$READ( 5, INFT ) (B 276) 

NC=4*NFO=I$JC=5U *PI =3ol^l3926326 $CNV=1.E-C 6 $ N M Y= 3 ( B 277) 

do; I S = : , 2*00118=1 ,5 *SH<1 )=«>3*3F ( 18) $SH(2 )=.7*BH< I B) $DGj IP=1,4 (B 278) 

002 I T=\ , S $RR = RP( IPI/BT ( IT) (R 279) 

S8 = SH( IS ) / 3 H < IB) $ Q S = i./S8$'*RITE(6?2)BF( IB) ,SB,RB ( B 280) 

FORMAT ( I H f ///hX4l-3/F = F7.3,10X4HS/B = F7.3,ltX*R/RETA = *F7.3,//5X3HY/ ( B 281) 
S , 7 X*DE L T A SLR 03*//) (B 282) 

D01IY=I f NMY$YH( IY, IS)=SH( IS)+FLOAT( lY-n/FLOATCNMY-l ) (B 28 3) 


IN THIS PROGRAM, THE INTEGRATION INTERVAL QNC IS A FUNCTION OF 
THE VARIABLE AYS ANALOGOUSLY TC ANC FCR THE SAME RE AS C NS THAT THE 
INTEGRATION INTERVAL QNC IS A FUNCTION CF THE VARIABLE ARGYS IN THE 
PRECEDING PROGRAM FOR DLTC3. 


AYS = YH( IY, I S)+SH< I S) $GNC=PI/ABS( AYS) *YS = YH( IY, IS) /SH( IS ) 
CALL INF INTO ( NQ, SMC, FG , NFQ , JQ , QNC ,CNV,QNT ) $Di = QNT 


(B 284) 
(B 2 8 5) 



APPENDIX B 


IF I Y = .! » YH=0 ANC AYS IS THE SAKE IN DOTH INTEGRALS. IF I Y =N MY » 
A YS=0 AND THE INTEGRAL IS ZERC. 


IF ( IY.EQ. 1 I GOTO 6 < B 286 ) 

IF I I Y. EC. AMY)G0TC4tAYS = SFUS)-YF( I Y , I S ) KNC= P I / ABS t A YS ) (B 2 87) 

CALL INF INTO! NO. SMQ. FC , NF Q , JQ , C NC ,C NV ,QNT > $D2 = QNT$GOT05 ( B 288) 

6 02 =01 $G0T05 (B 289) 

4 D2 = 0 » (B 290) 

5 DLTC3=-BS*<01+D2)/<2.*PI*PI)$WRITE<6,3)YS,DLT03 (B 291) 

PUNCH 7 , SB » BH ( IB) ,R3»YS*DL7G3 (B 292) 

7 FCRMAT(3<F7.3),Ffc«2»Fl('»6) (B 293) 

3 F0RMAT(4XF6.2,6XFI'i.6I (B 294) 

1 CGNTINUETSTCPJENC (B 295) 


SUBRCUTINE FUNCO IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST CVFR THE VARIABLE C IN ACCORDANCE WITH THE INTEGRATION 
SU8R0UT I NF MGAUSC. RfcFER TC CISCLSSION IN A PRECEDING PROGRAM. 

SUBROUTINE FUNC C < C , FC ) i C I MENS ICN FC (1) » S MR I 2 ) , FR ( 2 I , RNT< 2 ) , B F ( 5 ) (B 296) 
CONM0|\|SH ( 2 ),IS,YH(?,2),IY,RH,IE,SiPI,RB,AYS, I P $S = C $N R = 3 $ NFR = 2 ( B 29 7) 

IN THIS PROGRAM, THE INTEGRATION INTERVAL RNC IS MADE A FUNCTION 
CF THE VARIABLE RE ANC THE INCEXING INTEGER IP FOR THE SAME REASONS 
THAT THE INTEGRATION INTERVAL CNC WAS MADE A FUNCTION OF THE SAME 
VARIABLES IN THE PRECEDING PROGRAM FOR 0LT02. 

RNC-= C ',*RB / IPtCNV=i.F- A 7$JR=lCC 

CALLINFINTP(NR,SMR,FR,NFRfJR,RNC*CNV,RNT)$SYC=AYS*Q 
FO=SIN( SYC )* (0*R NT ( 1 ) +RNT (2)/<RB*RB*C))4RETURN$END 

SUBROUTINE F LNC R IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST OVER THE VARIABLES C AND R IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINES INFINTGt INFINTR, MGAUSC, AND MGAUSR. 


SUBROUTINE FUNC R < R , F R ) t D I ME N S I C N F R ( 2 ) , 8H ( 5 ) (B 301) 
COMMQNSHI ? ) , IS, YF< 3 ,2 ) , IY, BH, IB, Q, P 1 , RB, AYS (B 302) 
AR=SQRT(E*R/(R3*RR)+G*Q)$BR!=1./<RB*RB1+1. (B 303) 
FR(l)=i./< <BR)*R*R+C*C)*COSH( AR)*CCSMAR)-R*R) (B 304) 
FR (2 )=„./( ( BR.l+Q* 3/<R*R> ) *CC SH ( AR ) *CCSH< AR )-l. )t RETURN $END (B 305) 

1 1 NPT BH=.5,o75ti. *'. • 5,2, ,BT=1. ».8».fc».45,.3»RP=.l ,.45,2. ,7.5,$ (B 358) 


THIS PROGRAM MAKES USE OF TEE SUBROUTINE SUBPROGRAMS INFINTQ, INFINTR, 
MGAUSC, AND MGAUSR AS DEFINED AND DISCUSSED IN PRECEECING PROGRAMS. 


( B 298 ) 
18 299) 
(B 300) 
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FORTRAN PROGRAM FOR EVALUATING EQUATION (A-13) OF APPENDIX A. 

THE CALCULATIONS GIVE AN APPROXIMATION CLOSE TO THE SIMULTANEOUS 
ORIGIN OF THE THREE VARIABLES OF INTEGRATION TO TFE UPWASH INTER- 
FERENCE FACTORS FOR WINGS CENTER-MOUNTED IN PROGR ESS I VELY-MORE-OPEN- 
PERFOF'ATEO T t ST SbOTICNS. THE INTERFERENCE FACTORS CALCULATED 
REPRESENT TUG EFFECT OF HCRIZCNTAL BOUNDARIES ON THE INTERFERENCE 
VELOCITY POTENTIAL INSIDE TFE TEST SECTION DUE TO VERTICAL BOUNDARIES. 


IN ADUITICN TO THE CCMMCN VARIABLE NAMES ALREADY CEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TO THIS PARTICULAR PROGRAM. 

04 = THE UPwASH INTERFERENCE FACTOR CALCULATED. 

E P SI = A SMALL NUMBER, EPS ILCN-SUB-1 IN EQUATION (A-13). 

GPS2 = A SMALL NUMBER, E P S I L C N- SU 0- 2 IN EQUATION <A-i3). 

ALL OTHER VARIABLF N A v E S , EXCEPT AS NOTED, ARE USED MERELY AS 
BCCKKtFPlNG DEVICES. 

NAMELIST IN^LTS ARE BE, RT, RP 


PROGRAMWT ( INPUT , OUTPUT , TApF5= IN PL T,TAPC6 = OUTPUT, TAPE 1,TAFE2, PUNCH) 
DIMENSION *3H ( 5) » rfT ( 5 ) f RP(4 ) 

NAMELIST/INPT/MI ,BT ,F P$>READ( 5 , INFT) 

PI := Bol4iL c 2b3?^SZPSi = EPS' , = oJl$WR ITF (6,3) $wRITE(2,3) 

3 FORMUdH ,//5X3ha/H f 5X6HR/riE7A8X2hDA,//) 

DO ■ IB=;. f r <f DO; IP = I ,4$CC1IT = 3 , 5tRB = RP( I P ) / BT ( IT ) $RBH=RB/BU( IB ) 
RBS=R3U*F R UAPI = S QP T ( r PS2*E°S2+RBS ) 

AK = EPS2*AKL4Ri?S*ALGG(FFS2 + A«l)-ReS*ALrG(RBF)-HPS2*EPS2 
DA=FPS:*Rrt*AR/(2.*°I*PI*«RS)SV\RITE(i3f2)BH(ie),R9,C4 
w R I T j= ( 2 , 4 . ) R H { IB) ,RB,C4 
? FUr,viAT( 2( , XF7.3 , 2X) ,FE2o 3) 

PUNCF4, PH( I F ) , R B ,Q4 

4 FORMA T( 2( F 7. 3 ) ,Fu o 8) 

' CUNT I NUt*S TTPt^Nn 


$INPT BH=.3,o7^,2< 


, ,BT- V 


8 , o 6 f o 4 i 


R,RP = .l , o 4 5 , 2. ,7.5,* 


( 8 359 ) 
(B 360) 
{ B 361) 
<B 362) 
( P 363 ) 
IB 364) 
{ B 365 ) 
(B 366) 
IB 367) 
(B 368) 
(B 369) 
(B 370) 
( P 371 ) 
( B 372 ) 

(B 373 ) 
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FORTRAN PROGRAM FOR EVALUATING THE FIRST TRI P LE- I NTEGRAL PORTION 
OF EgiJATICN ( A- i A ) OF APPENDIX A. THIS CALCULATION BEGINS AT THE 
POINT THE IMMEDIATELY PPECFCING PROGRAM LEFT-OFF AND CARRIES THE 
INTEGPATICN THROUGH PART OF TEE REMAINDER OF THE INTEGRATION RANGE. 
THE CALCULATIONS GIVE AN APPROXIMATION AWAY FROM THE SIMULTANEOUS 
ORIGIN OF THE THREE INTEGRATION VARIABLES TO THE UPWASH INTERFERENCE 
FACT TP S EOF g ING S CENT ER-MCLNTFC IN PROGRE SSI VE L Y-MCRE-OPSN PER- 
FGRATLn TEST SECTIONS. TEE INTERFERENCE FACTORS CALCULATED REPRESENT 
THE EFFECT OF HORIZONTAL RCLNCARIES ON THE INTERFERENCE VELOCITY 
POTENTIAL INSICE THE TEST SECTION CUE TC VERTICAL BOUNDARIES. 


IN ORDER TC COMPARE WITH FCUATICN (A-14I OF APPENCIX A, NOTE 
THE FOLLOWING EQUIVALENCES. 

* 0 * OF THIS PROGRAM = GREEK LETTER *THETA* OF EQUATION (A-14). 
*P* OF THIS PROGRAM = GREEK LETTER *RHC* OF EQUATION IA-HI. 


IN ATCnifN TC THE COMMON VARIAPLE NAMES ALREADY DEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE 1C THIS PARTICULAR PROGRAM. 

OLT'K = T lJ F UFWASH INTERFERENCE FACTOR RESULTING ERC“ THE CAL- 
CULATIONS CF THIS PROGRAM, 

ALL OTHER VARIABLE NA W PS, EXCEPT AS MUTED, ARE USED MERELY AS 
BOCKKFf RING DEVICES. 


NAMELIST 1 N u L T S ARE BE, BT, RF 


PR OGK AM wT ( INPUT , CUT Ft T .TAPE 5=1 NPLT ,TAPE6=CLTPUT .TAPE 1 , TAP E2, PUNCH) (B 374) 
COMMUNISM <2 ) , IS, YF (3 ,: ) , I Y , BHI 5 ) , IE ,Q,R ,P I ,RB, AYS (B 3 75) 

DI MtNSIONF Q( ) , SMC ( : ) ,BT ( 5 ) ,RP (4 ) (B 376) 

NAMELIST/ INPT/BH,P.T ,RP$REAC(5, 1NPT ) ( 8 377) 

NFQ= :t u I = 3. 141 f.R 2<r 536 ENMI y = 3* AQ = C. $ BQ = P I /2. (B 378) 

DO!' IS = 1 ,2 TOOi. IR=: , 5 $SH(i ) = . 3*PP( IB) tSH<2 ) = . 7*BH< IB) *001 1 P = i ,4 { R 379) 

DO) IT = 1 , t <.RB=RP( IP) /HT( I T) $NC = 34 IFIXIPB/4. ) (B 380) 

SB = SH( IS ) /RF( 18 ) $BS=1. /S8?WR I TE ( 6 ,2 )Bh( I 8 ) , SB ,K3 (R 3811 

WRI TF(2,F ) 8H ( 18) ,SB,R8PC=RB*BS/(FI*PI*P1 ) (B 382) 

2 FORMAT ( 3 H .///• 3 X4P3/P = F7.2,1GXAHS/B = F7.3,ICX*R/BETA = *F7.3,//5X3HY/ (B 383 ) 

IS, 5X*0F LT A SUP 04*//) (B 384) 

DO] IY=l»NMYtYM IY,I S)=SH( IS)*FLOAT( I Y — 1 ) /FLOAT! NMY-1 ) ( B 385 ) 

YS = YH ( I Y , I S) /SH ( IS ) $C ALLMCAIJSG ( AC, BQ , NC, SMC , FQ.NFQ I $DL TO A = C * SMQ ( 1) ( B 3 8 6) 

wR I TE < 6, 3 ) YS,DL TC-i WRI TF ( 2 ,3 ) YS ,DLTG4 (B 387) 

3 FORM AT ( 2XF6.2 ,F XF ■ u .8 ) (B 388) 

PUNCH4 , SB » 8H ( I 3 ) ,RB,YS,DLT04 (B 389) 

4 FORMAT (3 (F7.“ ), Ffc.2 » F J f • • fc ) (B 390) 

i CONTI NUCiSTOPftND <B 391) 
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SUBROUTINE FUNCQ IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST OVER THE VARIABLE C IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINE MGAUSC. REFER TC DISCUSSION IN A PRECEDING PROGRAM. 


SUBROUTINE FUNC Q( « ,FQ ) $D I ME NS ICNFG (1) ,F R (1) , RNT (1) (B 392) 

CUMMCNSH<2),IS,YH(3,2),IY,BH<5), 1 8 ,S , R , P I , R 8 , AYS SS = Q SNR= I SNF R= 1 (B 3 93) 

JR = 5 ) SCNV=1 • F— ’ 5 SA Y S= YH ( I Y » I S ) +BH ( IB) $RNC = PI /ABS( AYS) ( B 394) 

CALLINFINTR(NR,FP,NFR, JR , R NC , CN V . RN T ) $C41=RN T ( 1) (B 395) 

IF< IY.EQ.J ) GO T06 $ A Y S= BH( I 8 ) - YH { I Y , I S ) S RNC = P I / ABS ( AY S > (B 396) 

CALLINFINTP(NR,FR,NFP,JR,RNC,CNV,RNT)ID42=RNT(1)$G0T05 (B 397) 

fa 042=041 (B 398) 

5 FQ (1 ) =CQS ( G ) * (D4 1 + C42 ) SRETURNJENC (B 399) 


SUBROUTINE FUNCR IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST OVER THE VAR I ABLtS C AND R IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINES INFINTR, MGAUSQ , AND MG AU SR . 


SUBROUTINE FUNC R ( R , FR ) $D I ME NS ICN FP ( I ) ,F P ( 1 ) , PNT ( I ) IB 400) 
CO^MONSHI 2 ) , I S,YH( 3 ,2 ) ,1 Y ,6H( 5 ) , IE ,C»T,PI , RB, AYS (B 401) 
T=R$NFP = : $JP = c ' SCNV=5oE-C5SPNC=PI/(2.*CnS<C) ) (B 402) 
IF ( PNC.GT. ( 2 . *P I ) ) PNC=2. *P I SNF = I F IX ( FNC*2./PI ) SIF (NP.LT. 1 )NP = 1 (B 403) 
CALLINFINTP(NP,FP,NFP,JP,PNC,CNV,PNT)$FR(1 ) =PNT( 1 ) $ RETURN SEND <8 404) 


SUBROUTINE FLNCP IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST CVtR THE VARIABLES F, C, AND R IN ACCORDANCE WITH THE 
INTEGRATION SUBROUTINES I \F I NT F » INFINTR, MGAUSP, MGAUSQ, AND MGAUSR. 


SLBRCLTINt FUNC F ( P ,FP) SCI BENS ICNFP (1 ) (B 405) 

CCMMC\SH(2 ) , IS, YM 3 ,2 ) , I Y, BH< S ) , IB ,Q,R,P I ,RB , AYS (B 406) 

RP=flH< IR > *P$SP=SH( I S> *PS A1 = SQRTI P^P^S INI Q)^S IN<Q)+R=^R) $A2=P*C0S (Q ) ( B 40 7) 

AR= 4YS*R$SNAR=AR$ IF (AR.GT.o05) SNAR = SIMAR) $CSAR=CCS( AR ) ( B 408) 

SHBP = BPS1F(BP.GT.o )5)SFRP = SINMBF)$CHEP=CCSH(BP) (B 409) 

AN=(P*SHBP*CSAR+R*C)'BP*SNAR)/(P*P+R*P)SSNG = Q (B 410) 

IFIC.GT.o • 5 ) SNO=SI N ( C ) SD? =SHP P * S F F F + R e *R B*SNC *SNQ *CHBP*C FBP (B 411) 

S HA. 1 = 1. $1 F ( A i o G T c o i 5 ) SHA 3 = S I NH ( A 1 ) / A } SCH A 1 = CC SH ( A 1 ) (B 412) 

D2=I»<8*P*SNC*SHA2)**2.+CHAl**2.tTRJ=SINH<SP)*AN/(EXP(8P)*Dl*02) (B 413) 

SNA2 = A2SIF(A t ;.GT.o • 5 ) SNA 2 = S I N ( A2 ) S C SA^ =C CS ( A2 ) $CSQ=CCS ( C ) (B 414) 

A3=( (RB*SNC)**L. >*P*SNA2*SHAT+CSC*CSA2*CHAJ (B 415) 

A4 = SHRP-k8*RR*SNC*SN<;*CHBP$A5 = A2*CSA2*SHA) -SNA2*CHA1 (B 416) 

FPU )=TR: *(EXP(BP)*A3+A4*A5) SRETURNSEND (B 417) 
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THE SUBROUTINE SUBPROGRAMS INFINTP AND MGAUSP IN THIS PROGRAM 
ARE NECESSARY TO ALLOW INTEGRATIONS OVER THREE VARIABLES RATHER THAN 
THE TWO VARIABLES CF PRECECING PROGRAMS. EXCEPTING FOR THE SUBROUTINE 
NAMES ANC WHATEVER VARIABLE NAME-CHANGES ARE NECESSARY TO ALLOW THE 
CALCULATIONS TO PRCCEEC WITHCLT CONFUSION, THEY ARE ICENTICAL WITH 
THF SUBROUTINES INFINTQ ANC MGAUSO AS DISCUSSED IN A PRECEDING 
PROGRAM. ALSO NOTE THAT TEE INITIAL VALUE CF THE LOWER LIMIT *AQ* 

HAS BEEN CHANGED FROM ZERO TO C.CI IN ACCORDANCE WITH THE LIMITS 
ON THE INTEGRALS BEING PROGRAMMED. 


SLBRCUT I Nb 1 NF I NT P ( NQ , FO F P , NF C , JQ , Q I NC , CCNV , Q INT ) 

DIMENSION SUMQM ) , FCFP(1 ) ,QINT(1) $0011=1, NFQ 

1 QI\T( I I = Jo $ AG=. 0 1$DC2 I 0=1 ,JQ$IF( IG-1 14, A, 3 
3 AQ=6U $GOTO Q 

A 8 Q = ) & 

9 + CSCALLMCA JSP( AO ,BQ,NG,SLMG,ECFP,NFC) 

CCCNV=.'. , NF G $CCQNV=GCCNV+ AGS ( SUMQ ( I ) ) 

5 QINT< I)=0 INT ( T ) + SUMC( I ) $ IF ( QCCN V . L7. C CNV ) GCTG6 $ I F ( IQ o E0. JG)GCT07 

GCTO? 

7 WRITE (6*8) SWRITE <? »6) SGGT06 

8 FORMAT ( 1H ,2X*!NTEGRAL NCNCCNV E R C- E NT VvITH GIVEN UPPER UNIT CN P*) 

2 CONTINUE 

6 RETURNEE NO 


SC8 ROUTINE MGAU SP < A ,E ,N ,SUM f FCFX t NUMBER ) 

DIMENSION U(5),R<5) ,SUMU) ,FOFX( 1 ) $00 ILL =i * NUMBER 

1 SUMLLl=U'’i 

IF (A. EQ.B) RETURN $U<I> = • 4 2 556 2 6 20 5 391 84 $U ( 2 ) =• 2 83 3H2 3029 8 h 376 

U ( 3 ) = © JL 6' 29 52l5 5 5 34e8$U(4>=.C6746£?16 6 55 5Ce$U5)=o0l30467 35741414 

R ( ll=. 147762) i2357376$R<2)=. 1 346 333 5 9654998 $R < 3 ) =„ 1 0 9543 1 81 2 5799 1 

R <4 > = .'>74725*6745 75 2 90 *R< 5 ) = . 03 3 3 3 5672 1 54 34 4 $F I NE=N 

OFLTA=F INE / < B*A) $003K=±,N$XI=K-1$FINF=A+XI / DEL T ASDG2 11 = 1,5 

UU = U( I D/OFLTA + F INEICALLFUNCPIUU ,FOFX ) $DQ 2 JO YBOY = 1 , NUMB E R 

2 SUM( J0Y30V) =RU I )*FCFX( JOY BOY) + SUM JCY8CY) $C03JJ = 1, 5 

UU= { 1. . -U( JJ > )/DELT A+FINE$CAI LFUNCP(UU,FCFX ) $D03NN = I , NUMBER 

3 SUM(NN) = R ( JJ ) *FCFX ( NNJ + SUMNN ) $CC7I JK = 1, NUMBER 
7 SUM( I J K ) = S U M ( IJK )/ DELTA $RE TUR N $ E N 0 


$ I NP T BH=.5».75»L. ,1.5,2 


BT-1# ,. P,#6,*A5,o3,PP-#I ,#45,2* ,7.5,$ 


{B 418) 
( B 419) 
( B 420) 
( B 421 ) 
(ft 422) 
(B 423) 
(B 424) 
(B 425) 
(B 426) 
(ft 427) 
(B 428) 
I B 429) 
(B 430) 


( B 443) 
(B 444) 
(8 445) 
( B 446 ) 
(ft 447) 
(B 448) 
(B 449) 
(B 450) 
(B 451) 
(B 452) 
( B 453 ) 
(B 454) 
(B 455) 

( B 482 ) 


THIS PROGRAM MAKES USE OF TFE SUBROUTINE SUBPROGRAMS INEINTR, 

MGAUSQ , ANn MGAUSK AS DEFINED AND DISCUSSED IN PRECEE C ING PROGRAMS. 
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FORTRAN PROGRAM FOR EVALUATING THE SECOND TR I PL E- I NTE GR AL PORTION 
OF EQUATION (A-14) OF APPENDIX A* THIS CALCULATION BEGINS AT THE 
POINT THE IMMEDIATELY PRECEDING PROGRAM LEFT-OFF AND CARRIES THE 
INTEGRATION ThROIGH THE REMINDER OF THE INTEGRATION RANGE# 

EXCEPT FOR TFT LIMITS OF INTEGRATION, THIS PROGRAM IS IDENTICAL 
WITH THE PROGRAM IMMEDIATELY PRECEDING. ONLY THE PORTION SHOWING 
THE CHANGE IN THE LIMITS OF INTEGRATION IS REPRODUCED HERE. 


SUBROUTINE FUNC R < R , F R ) $ 0 I ME N S IC NF R ( I > , F P ( 1 ) (B 509) 

COMMON SH (2 ) , IS, YH( 3 ,2 ) , IY,BH{ 5 ) , IE ,Q, T, PI ,RB, AYS ( B 510) 

SINCE T iE INTEGRATION LIMITS ON THIS INTEGRAL ARE FINITE RATHER 
THAN INFINITE AS KITH THE INTEGRAL OF THE PROGRAM IMMEDIATELY 
PRECEDING, THE SUBROUTINE INFINTP IS NOT USED IN THIS PROGRAM. 

T = R t.N F P = : $NP=;_t A P= * . $ D P= . i “1 $C ALL MGAU S P ( A P , 8 P , NP , FR ( 1 ),FP,NFP) (B 511) 

RCTUKN$END ( B 512) 

IN ACCORDANCE wITH THE INTEGRATION LIMITS OF EQUATION (A-lA) 

OF APPENDIX A, THE INITIAL VALUO FOR THE LOWER LIMIT AQ IN THE 
SUBROUTINE INFINTR IS SET AT C.C1 RATHER THAN ZERO. 


APPENDIX B 


FORTRAN PROGRAM FOR COLLATING DATA. THE PUNCHED-CARD OUTPUT 

OF the three preceding programs is the input tc this fpcgram. 

THE DATA ARE COLLATED, PRIMEC OUT AS A CHECK, THEN PUNCHED CNTO 
DA TA-PROCE SSI NG CARDS FOR LATER CCLLATIfN WITH TIE DATA REPRESENTING 
WINGS CF NTEW-RCUNTLC IN CLCSEC, RECTANGULAR, PERFORATED TE ST- SEC T I QNS 
AND WITH THE DATA REPKESNTING VARICLS COMPONENTS OF UFWASF 
INTERFERENCE FACTORS FOR WINGS C ENTER -MC LNTED IN PROGRESS I VEL Y- 
MCRP-OPE N , RECTANGULAR, PERFORATED TE ST -S ECT I CN S. THE COLLATED 
DATA LF THIS PPOOFAM GIVE TEE TOTAL COMPONENT FOR PROGRESSIVELY- 
MORE-OPfcN TFST-SECTIONS UF UPWASH INTERFERENCE FACTORS DUE TO THE 
EFFECT OF HORIZONTAL SCUNCA P I F S ON THE INTERFERENCE VELOCITY 
POTENTIAL INSIDE THE TEST-SECTION CUE TO VERTICAL BOUNDARIES. 


IN ADOITICN TC THE COMMON VARIABLE NAMES ALREADY CEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TO THIS PARTICULAR PROGRAM. 

04.1 = 04 CF TFE PROGRAM EVALUATING EQUATION (A-I3) OF APPENDIX A. 

D4 ? = DLT04 OF TFE PROGRAM EVALUATING TFE FIRST TR I PLE- 1 NTEGR AL 
PQPTICJN CF l- Q L A T I C N (A-J.4I CF APPENCIX A. 

043 = 0LTC4 OF TFE PROGRAM EVALUATING TFE SECOND TR I P L E- I NT EG R AL 
PORT ILN 0 F ECLAT ION (A-14) CF APPENDIX A. 

D43 = TFE RESULT OF COLLATION CF TFE FACTORS 041, D42 , AND D4 3 INTO 
THE UPWASH INTERFERENCE FACTOR REFFFSFNTING WINGS CENTER-MOUNTED 
IN PRDGRESSI VtlY-MCRE-OPEN, RECTANGULAR, PERFORATED T E ST- SEC T I ONS 
DLE TO TFF E FF C C T CF hCRIZCNTAL BOUNDARIES ON THE INTERFERENCE 
VELOCITY PUT ENT I AL INSIDE THE TFST-SECTICN OUF TO VERTICAL BOUNDARIES. 

ALL OTHER VARIABLE NAMES, EXCEPT AS NOTED, ARE USED MERELY AS 
BOCKKEEFING DEVICES. 


PRCGPAMViT ( INPUT , CUT FUT »T A F F 5 = INPUT , T A P F6 = CU T P UT , PUNC F. ) 

DIMENSION RH1 C T ') , BH2I6 Ul ,BF3 UCO I ,RB1 (1 GC I , RB2 (600 I ,R63 ( 6 )0 ) 
DIMENSION S P 2 ( S .G ) , SB ) (60 i ) * YS 2 ( 6PC ),YS3(t DC), 041(100) ,042(630) 
DIMENSION EU 3 (6 u. ) , C* 4 ( 6<' C ) , 04 5 ( 6 C L ) 

RE AC CSl ) ( BH1 ( I ) ,R el ( I ) , 041 ( I I , I = J , ir E ) 

1 FORMAT! 2(F7. i) ,FI'J, El 

R F AC ( 3 , 2 ) ( S (32 < I ) , HHZ < I ),;<«?( I >,Y$?< I) ,D4 2 ( I) , I =1 ,6'lC ) 

RF. ADIS,?) ( SB3(I ) ,PH3 ( I) ,RB3 ( I ) , Y S3 ( I ) , 04 3 ( I ) , I =1 , 600 ) 

2 FORMAT(3<F7.3),F6.2,FK.fc) 

wRITF ( 6 , " ) tWRlTE (6,i ) ( BH1 < I),REIIII,D41( 11,1 = 1,10') 

3 FORMAT (2X*B/ H*3X*R/RETA*3X MO LT41*//) 

WRI TE (6,4 ) $ WRITE (6 ,2) < SB2 ( 1 1 , BH2 (I) ,R82( I ) , YS2 (1) ,D42( I ) , 1 = 1 ,60C ) 
WRITE (6,8) tWRITE <6 , 2) < SB3< I I ,RH3 < I > ,RR3( II , YS3< I I ,043! II ,1 = 1 ,600 I 

4 PnRMAT(///2X«S/BM4X«e/H*3X*R/RETA*2XMY/SM4X*CLT42M// ) 
t> FORMAT ( / / /2XMS/BM4XMB/H*3XMR/BET A * 2 X M Y /SM4XMDLT43*/ / ) 

J = .itDCo IS = i , 2 $006 1 = 1,1 00 $CC6K=1, 3 t J = J + 1 $D4 4 ( J ) =D4 1 ( I ) 

6 D45< J )=044< J ) +042! J I+0F3IJI 

wRITE(6,7)$WKlTc(6f8)(SB2(IJ,SF3(I),RP3(I),YS3( 11,042(11 ,D43 ( I ) ,D4 
14(1 ),D45<I ) , I=). ,60.' ) 

7 F0RMAT(///2X«S/B*4X«B/II*3X*R/BETA*2X*Y/SM4X*QL T42 *5 X*DL T4 3* 8 X*DL T 4 
1 4* 5 X M OL T 4 5 * / /) 

8 rCRMAT(3(F7.3),F6,2,4(F10.6) ) 

PUNCH 2, ( SR3 ( I) , GH3 ( I I , R83 ( I I , Y S 2 (I) , 04 5 ( I) , I =1 , 6C J I 
STOP 1 END 
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APPENDIX B 


FORTRAN PROGRAM FOR EVALUATING EQUATION ( A-l ) ) OF APPENDIX A. 

THE RESULTS GIVE UPWASH INTERFERENCE FACTORS FOR WINGS CENTER-MOUNTED 
IN PROGRESSI VELY-MCRE-CPEN, RECTANGULAR, PERFORATED TEST-SECTIONS 
DUE TO THE EFFECT OF HORIZONTAL BOUNDARIES ON THE INTERFERENCE 
VELOCITY POTENTIAL OUTSIDE TEE TEST-SECTION DUE TO HORIZONTAL 
BOUNDAR IES. 


IN ADDITION TC THE COMMON VARIABLE NAMES ALREADY CEFINEO, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TO THIS PARTICULAR PROGRAM. 

0LT05 = THE UPWASH INTERFERENCE FACTOR CALCULATED 

ALL OTHER VARIABLE NAMES, EXCEPT AS NOTEC, ARE USED MERELY AS 
BOOKKEEPING DEVICES. 

NAMELIST INPLTS ARE RF, BT, RP 


PPCGRAMWT ( INPUT, CUT PUT ,T APE5 = IN PL T , TAPES = OU TP UT , TAP E I , TAPE 2 , PUNCH ) (B 578) 

COMMON SHI 2) , I $ , YH ( 3,2) ,IY,BH,I8,F,C,PI»Rfi,AYS,tP,BP,BPC,PQ,eQ,NMY (B 579) 

DIMENSIONFQ (1 ) ,SMQ( } ) f QNT (1 ) , EF( 5 ) , BT ( 5) ,RP (4 ) ( R 5 80) 

NAMELI ST/INPT/8H ,BT , R PSRE AD ( 5 , INPT) (B 581) 

NC = 1 $NFQ=jl tJC=5 ) $P I=3cI4j.5926E26 5CNV=2oE-C 2 $NMY=3t«NC=P I ( B 582 ) 

Dr)lIS = 2,2*DOiIH = 3,3$SH(J )=.3*EMie)$SM2) = .7*BH< I B ) $ DC2 I P = 3 , 4 (B 583) 

DO L I T = 1 , E $ R B - R P ( I P ) / B T ( I T ) (R 584) 

SB = SH ( IS) /6H{ IB) $RS=lo/S8$ WRITE (6 ,2 ) BF ( I P) ,SE,RB ( B 58 5 ) 

*RITE(2»Z)BHUBI ,SB»Re*C = eS/(PI*PI*PI) (B 586) 

2 FORMAT ( 1 H ,///oX4HB/h=F7.3,l 7 X4HS / B= F7 . 3 , i UX*R/ BET A= *F 7. 3 , //5X3HY/ (8 587) 

3S»5X*OELTA SUR 05*//) (8 588) 

00. U Y=: ,NMY$YH( I Y , I s ) =SH( IS ) *FLC AT ( IY-1 ) / FLOAT (NMY-1 ) ( B 589) 

YS = YH { I Y , I S ) / SH ( IS7 (B 590 ) 

CALL INF INTC(\C,SMQ, FQ,NFC, JO , C N C , C N V , QNT ) * CL T05=C*QN T< I ) (B 591) 

WP ITF ( 5, 1' ) YS,DL TC5$ WRI TE (2 ,3) YS ,CLTC5 <R 592) 

3 F0RMAT(2XF6o2,5XF 1 ,6) (R 593) 

pUNCHh,SE,BH(IB) ,RB,YS,DLTC5 {B 594) 

4 FORMAT {"{ F7 0 2 ), FE o2 ,F1C . E ) (B 595) 

.1. CONTINUt *5 TOP*; NO (8 596) 
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APPENDIX B 


SUBROUTINE FUNCQ IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTEREST OVER THE VARIABLF C IN ACCORDANCE WITH THE INTEGRATION 
SUBROUTINE' MGAUSCo REFER TC CISCUSSION IN A PRECEDING PROGRAM. 


SUBROUTINE FUNCC(Q,FQ1 $ C I MEN S IONFQ (11 , SMPIII ,FP( i ) ,PNT< 1) ,BH< 5) IB 597) 
COMMON SH ( 2 ) , I S, YH(3 ,2 ) , I V , BH , IB, F ,S, PI ,RB, AYS, IP, BP , BPQ, PQ, BQ,NMY ( R 59 8 1 
NP=5tNFP=! $PNC=5.*RB*FLOAT (I S ) /FLC AT ( I P ) SCN V=i . E-0 3 t JP=5C $ S=Q ( B 599) 
CALL INF INTP ( NP, S MP, FF ,NFP, JP ,FNC ,CNV,PNT ) $FQ ( i 1= G*PNT ( 1 ) (B 600) 
RETURNtEND ( B 601 ) 


SUBROUTINE FLNCP IS ONE WRITTEN TO EVALUATE THE INTEGRANC OF 
INTEREST OVER TFE VARIABLES P ANC C IN ACCORDANCE WITF THE INTEGRATION 
SUBROUTINES INFINTP, INFINTC, MGAUSP, AND MGAUSQ. 


SUBROUTINE FUNC P ( P , FP ) $D I ME NS ICNF P ID , F R ID ,SMR ( 1 ) , RNT (1) , BH ( 5 ) (B 602) 

COMMON SH< 2 ) , IS,YF ( j ,2 ) , IY, BH, IB, T,C,P I ,RB, AYS, IP, BP , BPQ , PQ , BQ , NM Y (B 603) 

T=P$NR=1 fNFP=3$JR=50$CNV=5.E-C4$eFQ=SQRT(P*P/(RB*RB)+G*Q) (B 604 ) 

BP=( io +lo / ( P B’i'R B ) ) *F*P$PQ=BP + C*Q1SQ=SH( I S ) * BPQ $BQ = BH ( IR ) *BPQ ( B 605) 

Ai. = S I NH{ SO) / trXP ( e J ) * EPO^BPQ ) (B 606) 

4YS=YHl IV, IS)+BH(IB)$PNC=PI/ABS(AYS) <B 607) 

CALL I NO INTRI NR, SMR , FP , NFR , JR ,R NC ,CNV , RNT ) $C5i =RNT ( I ) (B 608) 

IF ( IY.EQ„: ) GO TO 6 t AY S= BH ( IB l-YH IY , IS ) t RNC = P I / AB S ( A Y S ) ( B 609) 

CALL INF INTRINR, SYR , FR ,NFR , JR ,PNC ,CNV, RNT ) $C52=RNT (1 ) ( 6 63 0) 

GCTG5 ( B 611 ) 

6 D52 = 05’. (B 61.2) 

5 FP ( 1 ) = Ai * ( C5 i +11 5 2 ) $RdTURNtEND (B 613) 


SUBROUTINE FLNCR IS ONE WRITTEN TO EVALUATE THE INTEGRAND OF 
INTERTST CVLR TFE VARIABLES F, C, AND R IN ACCOR C ANC E WITH THE 
INTEGRATICN SLRRCLTINES INFINTP, INFINTC, INFINTR, MGAUSP, MGAUSQ , 
AND MG AU S R . 


S LBR CUT I N t FUNC R ( R ,FR ) tCIMF NS ICNFR ( 1 ) , BH ( 5 ) (B 614) 

COMMON SHI 2 ) ,1 S, YHI j ,2) ,IY,RH,IP,F,Q,RI,RB,AYS,IP,BP,BPO,FQ,BQ,NMY (B 615) 
ftPP = S(»RT ( P *P / I R 8 *R B ) + R*R ) t SR = SF ( IS ) *B P R $ B YR = A Y S *R $BP QR = B PQ*B PQ + R *R IB 616) 
PR=BP+R*ktA!=BPK*8PR/lBPCR*tPP*CCSHiePR)*CCSPiePR>-P*P}) (B 617) 

A5:=SIN(QI*C0SH(EPR)-C*C0S(Q)*SINP(BPR)/BPR (B 618) 

A5? = I3P3* ; C0S( 3YR) -RMSINIBYRItFRIl )=AX*A51*A52$RlTURN$£ND (B 619) 

& I NP T BH=.5,o75,lo , 1 .5,2o,BT=i. ».8,.6,.4&*.3,RF = .3. *.45,2. » 7 • 5 , t IB 695) 


THIS PROGRAM MAKES USE OF TFE SUBROUTINE SUBPROGRAMS INFINTP, INFINTQ, 
INFINTP, MGAUSP, Mf-AUSC, ANC MGALSP AS CEFINED ANC DISCUSSED IN 
PRECEDING PROGRAMS. 
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APPENDIX B 


FORTRAN PFCGFAM FCR COLLATING THE CAT A CALCULATED 8Y YEANS OF 
EQUATIONS ( A - 2 ) , (A-4), (A-7), ( A- 8 ) , (A-10), (A-13), AND (A-14) CF 

APPENDIX Ao THE PUNCHED-C A PD CUT PUT OF SEVERAL PRECECING PROGRAMS 
IS .THE INPUT TO THIS PROGRAY, THE DATA ARE COMBINED, PRINTED-OUT 
AS A CHECK, THEN PUNCHED CNTC C AT A- PROC E SS I NG CARDS HCR USE IN A 
PROGRAM WHICH GENERATES FINAL TABLES OF THE DATA, 


IN ACCITICN TC THE CCMMCN VARIABLE NAMES ALREADY DEFINED, 

THE FOl LOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM, 

DL T2 = THE UPWASH INTERFERENCE FACTOR DELTA-SUB-2 OF EQUATION (A-2) 
GF APPENDIX A* 

DLT3 = TEE UPWASH INTERFERENCE FACTCR DELTA-SUB-3 OF EQUATION (A-4) 
CF APPENDIX Ao 

DLT4 = TEE U^WASP INTERFERENCE FACTOR CELTA-SUB-4 OF EQUATION (A-14) 
OF APPENDIX Ao 

OLTf - THE UPWASH INTERFERENCE FACTCR CELTA-SUE-5 CF EQUATION (A-10) 
CF APPFKC ix A , 

DLTT = THE TOTAL UPWASH INTERFERENCE FACTCR, THE SUMMATION CF DLT2, 

D L T 3 , 0LT4, AND CL T6® 

ALL CTHFP VARIABLE NAMFS, EXCEPT AS NOTEC, ARE USED MERELY AS 
ROHKKELP ING DEV ICESo 


PROGRAMWTI INPUT, CL T PIT , PUNCH) 

DIMENSION SBC <3 ) ) , BHC (SO) ,YSC(3j),C2C( 30) ,D3C(3 )) ,OAC< 30 ) *D5C( 30) 
DIMENSION S B U ( 6 . ' , 4 ) ,BHO( fc(\) ,A ) , RBC (ECO, A } , YSCMoOJ, A) ,DLT (6C 0,4) 

DIMENSION Q L T 7 ( 6 < , ) , CL T? ( 6v. 0 ) , CL T A ( 6( < ) , DL T 5 ( 60U) ,DLTT( 6C0 ) 

READ v , ( SBC ( I > ,3HC( I ) ,YSC( I) ,D2C( I) ,D3C( I) , DAC ( I ) , 05 C ( I > , I = 1 , 30 > 

? F0RMAT(?(F 7 «3),F6oC,4(F?2oE ) ) 

PRINT? 

2 FORMAT (2X*S/3*4X*B/F*3X*Y/S*5X*DLTC2*7X*DLTC3*7X*DL TC4*7X*DLTC5*/) 
PR INTI , { SBC ( I) ,BHC( I ) , YSC( I I ,D2C ( I ) ,D3C( I ) , CAC( I ) ,D5C( t ) ,1=1 ,30) 
kC A C? , ( ($RO ( I f J ) ,BHC( I ,J ) ,RB0 ( I , J ) , YS0( I , J ) ,DLT< I , J ) , 1 =1,600 ) , J = I, 

i 4 ) 

3 FORMAT ( 3 ( F 7 , 3 ) , F6,2 , F 1 C ,6) 

0 0 A J = 1 , 4 $ L = J + 1 

D R I N T 5 , L , (SBO(I,J),eHn(I,J),ReC( I , J ) , Y SO ( I,J),DLT(I,J) , I=1,6CQ) 
b FL'RMAT{ ///2 x*S/B*4X*Q/H*2X*R/8ETA*2X*Y/S*?X*CELr A SUB 0*Ii f //(3(F7 
lo 3) , F 6 o 2 , F 1 C o 6 ) ) 

A CONTINUE 

J = $0061 = ? , I U$MC=( I-I )*3 + ! i>NC=MC + 2$DC6lRB=l ,2 » $006 I Y =MC , NC $ J = J+l 
DLT2( J)=n?C( IY)+CLT (J f i) $CLT3( J) = C3C( IY)+CLT(J,Z) 

0LT4( J )=|U.C ( I Y)+CLT( J , 3) $DLTS( J) =C5C( IY) + CLT< J ,4) 

6 DLTT<J)=DLT2(J)+CLT3<J)+CLT4(J)+CLT5(J) 

PR INT7,{ SBC( I ) ,B U C U ,1 ) ,RBO( I , YSC ( I ,1) ,CLT2 ( I ) ,DLT3 ( I ) , CLTA ( I 

1 ) , DLT 5 ( I ) , 0 L T T ( I ) , K. ,5u. ) 

7 FORMAT(///2X*S/il*4\*t3/H*2X*R/BETA*3X*Y/S*3**CELTA SUB 2*2X*DELTA S 
1UB 3*3X*OtLTA SUP a*:y*dELTA SLR 5*3X*TOTAL DE L T A* / / ( 3 ( F 1 o 3 ) , F 6, 2 , 

2 5 ( ?XF8o6,?X) ) ) 

PUNCH R , (SB°( I, . ) ,BHL( 1,3 ) , RBO ( 1,1 ) ,YSG( 1,1 ),DLT2( I ) ,DLT3( I ) ,DLT4(I 
1 ) ,OLTS( I ) , )LTT( I ) , I = I , F k >. ) 

8 FORMAT ( 3 F 7 o 7 , F & o 2 , 6 F ’ ■ o 6 ) 

S TCP $ t NO 
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APPENDIX B 


TCKTRAN PROGRAM FOR GENERATING TABLES CF THF PREVIOUSLY CAL- 
CULATED ANT COLLATED DATA, TEE DATA CLTPUT FROM THE IMMEDIATELY 
PRECEDING PROGRAM ARE THE INPLT FCR THIS PROGRAM* THIS PROGRAM 
ARRANGES T He DATA AND PRIMS TEEM IN TABULAR FORM* 


IN AOOITICN TC THE COMMON VARIABLE NAMES ALREADY DEFINED, 

THE FOLLOWING NAMES ARE OF IMPORTANCE TC THIS PARTICULAR PROGRAM. 


D2 

= D LT 2 

OF 

THb 

IMMEC IAT ELY 

PRECE C ING 

PROGRAM 

D 3 

= DLT2 

OF 

THf 

IMMEC I ATELY 

FRECECING 

PROGRAM 

04 

= OLTA 

OF 

TMt 

IMMCC I ATELY 

FRECfC ING 

PROGRAM 

D5 

= LLTb 

or 

THE 

IMMbCIATFLY 

PRECEDING 

PROGRAM 

DT 

= DLTT 

CF 

THE 

IMMEC I ATELY 

FRECECING 

PROGRAM 


ALL OTHER VARIABLE NAMES, EXCFFT AS NOTEC, ARE USED MERELY AS 
BOOKKt f PING DLVI Ct S. 


PROGRAM V.'TIN( INPLT, CLTPUT) 

DI MENS I ON SB <6 * ) ) , BHI60L ) , RB ( 6 30 ) , R < 6l ? ) , Y S ( 600) ,D2< 60C ) ,0 3( 600) 

D I MENS I Of' DM6' o ), 05(6^0) ,DT(6CC) , BET A (600 
DIMENSION X(2,,5,^),Y(?0,5,2),XP(22),YP<?2) 

Rb ad;. , ( SEMI ) ,BH( I ) , RP ( I ) , YS( I ) ,D 2 ( I ) ,03 ( I ) , DA (I) , D5 ( I ) , DT ( I ) ,1=1,6 
] ) S N B H = - 2 9 $ N P = - ] MNP ET=NYS = -2 

1 FORMAT ( 3 0 7* 2 , F 6 © 2 » 5 F i 0 * 6 ) 

D02 I S = 1 , ; $DC? IDH = \ , 5*NRB=~SDCMPM , 2 $N 6 H=N P E + T 0 
PR I N T 5 , B H ( NBH) , SBC NRH) $006 I R = ) ,1C$NBET=N9ET+3$NR8=NRR+1 
6 FOKMATH H*. /ftX^TUNNOL WIDTE-TC-LEICET RATIO, P /H= *F6 • 2 , P X * * I NG- SP AN 
l -TOTUNNEL-to IDTH RATIO, S /B = *F 6 • 2 / / /4X*R/BET A* 8X*Y/ S*<= X*C EL T A*8X*0 
2ELTA*BX*DElTAMX*ntL7 A*R X*TOT AL * 5 X* DEL TA /D EL T A * / 3 I X*SUB 2*8X*SUB 3 
3*8X*SUR ^**8X*5J8 5*FX*DELTA*RX*SLR C*//) 

PRINT?, RB(NBfcT) 

Y(NKB, IBH, I S ) =0T (NBET)$X(NRB» IEh,IS)=RR(N8ET) £MYS=NBET+2 
PKINTiOt ( YS ( I) ,02( I ) ,D3( I) ,DM I) ,05 (I) ,DT< 1) , I=NBET,MYS) 

6 CUNT I NUF 
A CONTINUE 

PRINT 11 ,( X( I yIRHfIS) ,Y(I , I BH , I S ) ,1=1, NRB) 

CONTINUE 
CONTINUE 

FORMAT ( I H + , 1 X F 7 •? ) 

1 j FORMAT! 1H ,l6XP4*),5(7XF6o3)/]7XF4©I,5(7XF6.3)/17XF4*l,5(7XF6.3)/) 
U FORMAT (1H3 , 1 X*k/ BET A*6X*TGT AL * / 1 4 X *DL L TA */ / < 2 XF 7* 3 , 3 XF 6. 2 ) ) 

STOP SEND 


(B 750) 
(B 751) 
(B 752 ) 
(B 753) 
(R 754) 
(B 755) 
(B 756) 
(B 757) 
(B 753) 
<B 759) 
(B 76 3 ) 
( B 761) 
(B 762) 
<B 763) 
(B 764) 
(B 765) 
(B 766) 
(B 767) 
<B 768) 
(B 769) 
(B 770) 
(R 771) 
(B 772) 
(B 773) 
(B 774) 
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TABLE I.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING- SPAN-TEST- SECTION- WIDTH RATIO s/b OF 0.3 

(a) Tunnel width-height ratio b/h of 0.50 


Approximated 6 


R /0 

y/s 

6 2 

S 3 

S 4 

6 5 

<?2 + 6 3 + 6 4 + S 5) 

CLOSED 

0.0 





.266 

.100 

0 .0 

.205 

.027 

-. 014 

-. 002 

.216 


. 5 

.208 

.027 

-.014 

-.002 

.219 


1 .0 

.216 

. 026 

-.014 

-. 002 

.227 

. 125 

0.0 

.194 

.026 

-. 012 

002 

.206 


. 5 

.197 

.326 

-.012 

-.03 1 

.209 


1 .0 

.205 

.025 

-.012 

00 1 

.217 

.167 

0 .0 

.179 

.023 

-.010 

001 

.191 


. 5 

.180 

.323 

-.010 

-.001 

.193 


1 .0 

.188 

.023 

-.009 

001 

.201 

.222 

0.0 

.159 

.020 

007 

-. 000 

.172 


. 5 

. 161 

.02 0 

-.007 

-.000 

.174 


1.0 

.169 

.020 

-.007 

. 000 

.181 

.333 

0.0 

.127 

.015 

005 

. 002 

.139 


. 5 

. 129 

.015 

-.005 

.002 

. 141 


1 .0 

.135 

.014 

004 

. 002 

.147 

.450 

0.0 

.103 

.009 

-.003 

. 003 

.109 


.5 

. 101 

.009 

-.003 

.003 

.110 


1.0 

. 107 

.009 

-.003 

. 003 

.116 

.56 2 

0.0 

.079 

.004 

-.002 

. 004 

.0 84 


. 5 

. 080 

.004 

-.302 

.004 

.086 


1.0 

.084 

.004 

-.002 

. 005 

.090 

.750 

0.0 

.049 

-.004 

-. 001 

-00 6 

.051 


. 5 

. C 50 

-.004 

-.001 

.006 

.052 


1 .0 

.053 

-.004 

-.001 

. 00 7 

.055 

1 .000 

0.0 

.021 

-.0 12 

.000 

.00 9 

.017 


. 5 

. 021 

012 

.300 

.009 

.018 


1 .0 

.023 

-.012 

. 000 

. 009 

.019 

1 .500 

0.0 

-.0 17 

-.024 

.004 

. 01 1 

026 


a 5 

-. C 1 7 

-.024 

.004 

.012 

- .026 


1 .0 

-.018 

- .024 

. 004 

. 012 

-.027 

2.000 

0.0 

-.041 

-.032 

.008 

. 01 3 

-. 052 


. 5 

-.041 

-.032 

. 008 

.013 

- .052 


1 .0 

-.044 

-.032 

.007 

. 013 

-. 054 

2.500 

0. 0 

-.056 

-.033 

.011 

.014 

- .069 


.5 

-.057 

-. 037 

. 01 1 

.01 4 

-.070 


1 .0 

-.061 

-.037 

.01 1 

.015 

073 

3.333 

0. 0 

-.074 

-.344 

.015 

.015 

- .087 


.5 

- . C 75 

-.044 

. 01 5 

.015 

-.0 8 8 


1.0 

-.079 

-.043 

.014 

.016 

-. 093 

4.444 

0.0 

-.087 

-.349 

.020 

.01 6 

- .100 


.5 

- . C 89 

-. 048 

.020 

. 016 

-.102 


1.0 

-.094 

-.048 

.019 

.017 

107 

6.667 

0. 0 

-.102 

-.354 

.024 

.017 

- .116 


.5 

-. 104 

-.054 

. 024 

. 01 7 

-.118 


1. 0 

-.111 

-.053 

.023 

.017 

- .124 

7.500 

0.0 

-. 106 

-.055 

.025 

.017 

-.119 


.5 

-.108 

-.055 

. 025 

. 017 

-.121 


1. 0 

-.114 

-.055 

.024 

.017 

- .127 

5. 375 

0. 0 

-.111 

-.057 

.027 

.017 

-.124 


.5 

-.113 

-.057 

. 027 

.017 

-.126 


1.0 

-.120 

-.056 

.026 

.017 

- .133 

l 2. 500 

0.0 

-. 117 

-.059 

.029 

.017 

-.130 


.5 

-.119 

-.059 

. 025 

. 01 7 

-.132 


1.0 

-.127 

-.058 

.028 

.018 

- .139 

1 6. 667 

0 . c 

- . 121 

061 

. 031 

.017 

-.134 


.5 

-.124 

-.060 

.031 

.017 

-.136 


1.0 

-.131 

-.060 

.030 

.018 

-.144 

25. 000 

0.0 

-.125 

062 

. 033 

.017 

-.138 


140 

-.148 
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TABLE I.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN-TEST-SECTION- WIDTH RATIO s/b OF 0.3 - Continued 

(b) Tunnel width-height ratio b/h of 0.75 





TABLE I.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN— TEST-SECTION- WIDTH RATIO s/b OF 0.3 - Continued 

(c) Tunnel width- height ratio b/h of 1.00 


Approximated 6 


R/0 

y/s 

6 2 

6 3 

6 4 

5 5 

( 6 2 + 6 3 + 6 4 + 6 J 

CLOSED 

0.0 





.138 

.100 

0.0 

.102 

.053 

-.042 

.004 

.117 


. 5 

.104 

.052 

-.042 

.004 

.118 


1 .0 

.103 

.049 

— • 04 C 

.004 

. 121 

.125 

0. 0 

.097 

.050 

-.039 

.004 

.112 


.5 

.098 

. 049 

-.03 9 

.004 

.113 


1 .0 

. 103 

.046 

- .0 3 B 

.004 

. 115 

.167 

o 

o 

.089 

.046 

-.036 

.004 

.103 


.5 

.090 

. 045 

-.03 5 

.004 

.104 


1.0 

.094 

.042 

-.034 

.00 5 

.107 

• 22 2 

0.0 

. C 79 

.040 

-.031 

.004 

.092 


.5 

.081 

.039 

-.031 

, 005 

.093 


1.0 

.084 

.036 

-.030 

.005 

.096 

.333 

0. 0 

. C 63 

.029 

-.025 

.005 

.072 


.5 

.064 

.028 

-.024 

.005 

.073 


1.0 

. 068 

.026 

-.024 

.006 

.075 

.450 

0. 0 

. C 50 

.017 

-.020 

.005 

.053 


.5 

.051 

.017 

-.019 

.006 

.054 


1.0 

.053 

.015 

-.019 

.006 

.056 

. 562 

0. 0 

. C 3 5 

.007 

-.015 

.006 

.037 


.5 

.043 

.007 

-.015 

.006 

.037 


1.0 

. 042 

.005 

-.015 

.007 

.039 

. 750 

0. 0 

. 025 

-. 008 

-.010 

.006 

.014 


.5 

.025 

-.008 

-.010 

.007 

.014 


1.0 

. 027 

-.009 

-.010 

.007 

.015 

1. 000 

c.o 

.010 

-.024 

-. 004 

.00 7 

-.011 


.5 

.011 

-.024 

-.004 

.007 

-.011 


1.0 

. 011 

-. 025 

-. 004 

.008 

- .0 10 

1.500 

0.0 

-.009 

-. 048 

.005 

. 00 8 

-.044 


.5 

-.009 

-.048 

.005 

.008 

-.044 


1.0 

- . CC 9 

-. 048 

.004 

.009 

-.043 

2.000 

0.0 

-.020 

-.064 

.011 

. 00 8 

-.065 


. 5 

021 

-.063 

.011 

.008 

- .065 


1.0 

-.022 

-.063 

.011 

. 009 

-.064 

2.500 

0.0 

-.023 

-.075 

.016 

.00 8 

-.079 


. 5 

-. 029 

-.074 

.016 

.008 

-.079 


1 .0 

-.030 

-.073 

.015 

. 00 9 

-.079 

3.333 

0.0 

-.037 

-.087 

.022 

.00 8 

-.094 


. 5 

C 37 

-.086 

.021 

.008 

- .094 


1 .0 

-.040 

-.084 

.021 

. 009 

-.094 

4.444 

0.0 

-.044 

-.096 

.028 

.00 8 

-. 1 C 5 


. 5 

-. C 45 

-. 096 

.028 

.038 

- .104 


1 .0 

- .047 

-.093 

.027 

. 009 

-.104 

6.667 

0.0 

-.051 

-.107 

.034 

.007 

-. 117 


.5 

- .052 

-.106 

. 033 

.007 

-.117 


1 .0 

-.055 

-.103 

.033 

.009 

-. 117 

7.500 

o 

o 

-.053 

-.109 

.035 

.007 

- .120 


.5 

- .054 

-.108 

. 03 5 

.007 

-.120 


1.0 

-.057 

-.105 

.034 

.009 

- .119 

5.3 75 

0.0 

-. 056 

-.113 

.038 

.007 

- .124 


.5 

- .057 

-.112 

.038 

.007 

-.124 


1.0 

-.06 0 

-.109 

.037 

.00 9 

-.123 

1 2.500 

0.0 

-. 058 

-.117 

.041 

.006 

- .128 


.5 

-.060 

-.116 

.041 

.00 7 

-.128 


1.0 

-.063 

-.112 

.040 

.00 8 

- .128 

1 6.667 

0.0 

-. 061 

-.120 

.043 

.006 

-.131 


.5 

-.062 

-.119 

.043 

. 007 

-.131 


1.0 

-.066 

-.115 

.042 

.008 

- .131 

2 5. 00 0 

C.O 

-. C 6 3 

-.123 

.045 

.006 

- .135 


• 5 

-.064 

-.122 

. 045 

.006 

-.135 


-.134 


136 





























TABLE II.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN-TEST-SECTION- WIDTH RATIO s/b OF 0.7 

(a) Tunnel width -height ratio b/h of 0.50 



TABLE H.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN- TEST-SECTION- WIDTH RATIO s/b OF 0.7 - Continued 

(b) Tunnel width -height ratio b/h of 0.75 








Approximated 5 

R //3 

y/s 

5 2 

6 3 

5 4 

S 5 

( Q 2 + 6 3 + 5 4 + 6 5 ) 

CLOSED 

0.0 





.191 

.100 

0.0 

.148 

.038 

-.028 

.001 

. 159 


.5 

. 164 

. 036 

-.027 

.001 

.173 


1 .0 

.230 

.0 30 

-.023 

. 001 

.238 

.125 

0.0 

.141 

.036 

-.026 

.001 

. 152 


.5 

.156 

. 034 

-.02 5 

.001 

.166 


1 .0 

.220 

■ .02 8 

-.021 

.002 

.228 

.167 

0.0 

.129 

.032 

-.023 

.002 

.141 


.5 

.144 

. 030 

-.022 

. 002 

.154 


1.0 

.204 

.025 

-.019 

.002 

.213 | 

.222 

0.0 

.116 

.328 

-.019 

.002 

.127 


.5 

.129 

.026 

-.oie 

. 002 

.139 


1.0 

. 186 

.022 

-.016 

.003 

.194 

. 333 

0.0 

.093 

.0 20 

-. 015 

.003 

-102 


.5 

.104 

.018 

-.014 

.004 

. 112 


1.0 

. 153 

.014 

-.012 

.004 

.160 

. 450 

0.0 

.074 

. 012 

-.011 

.004 

.0 78 


.5 

.083 

.010 

- .011 

.005 

. 068 


1 . 0 

. 124 

. 007 

-.009 

.006 

.128 

.562 

0.0 

.058 

.004 

-.009 

. 005 

.059 


. 5 

.066 

.003 

-.008 

.006 

.067 


1.0 

. 100 

. 000 

-. 0 C 7 

.007 

.100 

.750 

0 .0 

.037 

- .006 

-.006 

. 007 

.032 


. 5 

.04 2 

-.007 

-.005 

.00 7 

.037 


1.0 

.066 

-.009 

-.004 

.009 

.061 

1.000 

0.0 

.016 

-.018 

-.002 

. 00 8 

.004 


.5 

. 013 

-.019 

-.002 

.009 

.007 


1 .0 

.030 

-.020 

-. 001 

. 011 

.020 

1.500 

0.0 

-.012 

-.036 

.005 

.010 

-.033 


. 5 

-.014 

-.036 

.005 

.011 

-.034 


1.0 

-.019 

- .035 

• 0 C 4 

. 013 

-.036 

2.000 

0.0 

-.030 

-.047 

.011 

.011 

-.056 


. 5 

-.034 ! 

-.047 

.010 

.012 

-.058 


1 .0 

-.051 

- .045 

. 009 

. 01 5 

-.071 

2.500 

0.0 

-.042 1 

-.055 

.015 

.012 

-.070 


.5 

- .048 

-.054 

. 01 5 

.013 

-.074 


1 .0 

-.073 

-.051 

.013 

.016 

-.095 

3.333 

0.0 

-.055 

-.064 

.021 

.012 

- .086 


.5 

-.063 

062 

. 02 C 

.013 

-.092 


1 .0 

-.097 

-.059 

.018 

. 017 

121 

4.444 

0 . c 

-. 065 

-.071 

.028 

.012 

- .097 


.5 

-.075 

-. 069 

. 02 6 

. 014 

-.104 


1 . 0 

-.117 

-.065 

.023 

.018 

- .141 

6 . 667 

0 . 0 

-. 077 

-.078 

.034 

.012 

-.109 


.5 

-.008 

-.076 

.032 

.014 

-.118 


1 . 0 

-.133 

-.071 

; .028 

.019 

- .163 

7 . 500 

0 . 0 

-.079 

080 

.035 

.012 

-.112 


.5 

-.091 

-.078 

.034 

.014 

-. 122 


1.0 

143 

-.072 

.029 

.019 

-.168 

5 . 375 

0.0 

-.084 

-.083 

.038 

.012 

-.116 


.5 

-.096 

-.080 

.036 

.014 

-. 126 


l.C 

-. 151 

-.075 

.031 

.019 

-.175 

12. 500 

c.o 

-.008 

-. 085 

. 041 

• Oil 

-.121 


.5 

-.101 

-.083 

.039 

.014 

-.131 


1.0 

-.160 

-.077 

.034 

.02 0 

-.184 

16.667 

0.0 

-.091 

-. 087 

. 043 

.011 

-.124 


.5 

-.105 

-.085 

.041 

.013 

135 


1 . 0 

-.166 

079 

.035 

.020 

-.190 

25 . 00 C 

0.0 

- .095 

-. 090 

. 046 

.011 

| -.128 


.5 

-.109 

-.087 

.044 

.013 

-. 139 


1.0 

-.173 

-.080 

.037 

.020 

-.196 

OPEN 

0.0 





128 


53 




















TABLE n.- UPWASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN— TEST -SECTION- WIDTH RATIO s/b OF 0.7 - Continued 

(d) Tunnel width -height ratio b/h of 1.50 


Approximated 6 
( 6 2 + ®3 + S 4 + 6 5 



TABLE H.- UP WASH INTERFERENCE FACTORS 6 IN RECTANGULAR TEST SECTIONS 
CALCULATED BY APPROXIMATION PROCEDURE WITH VARYING PERMEABILITY 
FACTORS FOR WING-SPAN— TEST-SECTION-WIDTH RATIO s/b OF 0.7 - Concluded 

(e) Tunnel width -height ratio b/h of 2.00 














